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PREFACE 
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INTRODUCTION 


Since 1963, an extensive investigation of the inducer flows has 
been carried out at the Department of Aerospace Engineering at The 
Pennsylvania State University under NASA sponsorship. A model of an 
inducer was constructed and characteristics of two, three and four bladed 
inducers were investigated. The inducer blade had radially variable 
thickness, camber, and tangential extent of average 290 degrees. The 
flow though a long and narrow passage is subjected both to the potential 
flow effect determined by the rotor blade geometry and to the effect 
of the viscosity of fluid. Attempts to predict the flow in these 
passages is dealt with in Refs. 10 and 15. 

It is seen that, due to the effect of viscosity, the fluid in 
the immediate vicinity of the blade surfaces is set in rotating motion, 
and is subjected to centrifugal and Coriolis forces. As a result, 
direction of the flow outside the boundary layer becomes different from 
the direction of the flow inside the boundary layer, i.e. , the skewed, 
or three-dimensional, boundary layer is formed. Thus, it becomes 
essential to investigate the behavior of three-dimensional boundary 
layer flow in order to understand behavior of. the flow in the inducer 
passage. 

Since the blade passages are long and narrow, the boundary layers 
on the inducer surfaces experience channel effects. They interact 
strongly with the comparatively inviscid external flow, and with the 
annulus wall boundary layer, producing the observed radial inward flows. 
(23,24). These radial flows are found to increase as the blade boundary 



2 


layers grow and merge with stress-reversal, and reduction in mainstream 
velocity, etc. (23, 39). Since these phenomena have not been observed in 
a stationary channel, these are anticipated to have been caused by the 
rotation of the blade passages. 

General treatment of the three-dimensional boundary layer flow is 
extremely difficult because of the mathematical difficulties involved 
in solving the governing equations. The method used most often to solve 
this problem is to assume the mean velocity distribution in the boundary 
layer based on various laws that are obtained mainly for two-dimensional 
boundary layer flows. Though this approach has been successful for the 
treatment of the three-dimensional boundary layer flows that do not 
deviate much from the two-dimensional ones, it is expected that such 
approach cannot be applied satisfactorily to the highly skewed boundary 
layer flows which are created by the effect of rotation, i.e., the 
rotating type three-dimensional boundary layer flows. With an intention 
to obtain general information on the rotating type of boundary layer 
flow, a helical blade was constructed and some characteristics of 
turbulent boundary layer flow on the rotating helical blade were investi- 
gated (5,9). The previous analysis of the flow carried out at Penn State 
are based on assumptions for velocity profiles, which is largely derived 
from the existing information on two-dimensional and three-dimensional 
boundary layers. The prediction so obtained are only qualitative. The 
main reason for this being the lack of information on the velocity pro- 
files and shear, stress in a narrow rotating channel of type used in an 
Inducer. The objective of this investigation is to provide this important 
scientific information in the eventual viscid solution of the flow field 


in an inducer. 
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An attempt is made in this report to predict the three-dimensional 
boundary layer characteristics in a rotating channel. Sections 2.2 and 
2.3 deal with the velocity profile analysis and the momentum integral 
solution is described in section 2.4. An approximate solution of the 
flow in the radial flow reversal region (mid passage) is dealt with in 
section 2.4 and a future method of approach is outlined in section 2.5 
In order to investigate the effect of the fluid viscosity, a rotor 
assembly was constructed with a simple geometry. The fabrication details 
of this new rotor assembly, machined accurately to avoid vibration, is 
described in section 3.1. Instruments used are described in section 3.2. 
Measurements of the blade static pressure, limiting streamline angle and 
skin friction were carried out on the blade surfaces of open throttle 
conditions (Sec. 3.5, 3.6, 3.7). Total pressure and flow angle measurements 
were carried out inside the rotor channel (Section 3.8). A throttle 
was installed to study the flow in the inducer at different flow coefficient 
(Section 3,3) and the flow measurements at the rotor exit are reported 
in Section 3.4. All measurements were carried out at 450 RPM and at 
Reynolds Number (based on tip radius) 6.6 x 10^. Unless otherwise 
indicated, the measurements refer to no throttle conditions. 

2. THEORETICAL CONSIDERATIONS 

* Flow in the rotating helical channel can be divided into two parts. 

One is the boundary layer flow near the blade surfaces 

where the effect of viscosity of the fluid is large, and the other is 

the flow outside the boundary layer where the effect of viscosity of the 
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fluid is negligible. Since the boundary layer flow has a rotational 
motion, it is subjected to two additional forces, centrifugal and 
Coriolis forces. As a result, a radial flow, or cross flow, is pro- 
duced. The direction of the boundary layer flow becomes different 
from the direction of the flow outside the boundary layer, i.e., three- 
dimensional boundary layer is formed on the blade surfaces. 

The first general model of the three-dimensional boundary layer 
flow was proposed by Prandtl in 1946 (16). He considered that the flow 
in the free stream direction u and the cross flow w can be represented 
by 


u_ 

U 


c(n) 



(i) 


where g, G are universal functions of r) = z/6. U is a free stream velocity, 
is the limiting streamline parameter, z the distance normal to the 
surface, and 6 the boundary layer thickness. This model has been used 
by many. Perhaps the roost well known model was given by Mager (11). 

This model can be considered as an extension of two-dimensional power 
law approximation to three-dimensional boundary layer flow. 

Coles (1) suggested that his model for the two-dimensional boundary 
layer flow could be extended to the three-dimensional boundary layer flow. 

He considered that the mean velocity profile could be represented by. 


u = U T f(z + ) + U^nwCn) , where r) = z/6 . 


( 2 ) 
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where f, W are universal functions representing the "law of the wall" 
and the "law of the wake" respectively. U is a friction velocity 
vector; its magnitude is given by /t /p and its direction is the same 
as the wall shear stress direction. 7T is a tensor such that u is in 
the direction of the free stream velocity at z = 6. Validity of the 
Coles model was experimentally investigated by Homung and Joubert (4). 
Pierce (1J) also demonstrated that the Coles model was appropriate for 
the skewed flow. 

In 1960, Johnston (7) presented a new approach to the three- 
dimensional boundary layer flow analysis. He considered that the cross 
flow component was related to the main flow component by 


w 

U 


£ ( v 



(3) 


where e Q , A are parameters. In the outer part of the boundary layer 
he obtained a linear relationship given by 


If = A(l-J) (4) 

By considering that in the region very close to the wall the direction 
of the flow must be the same as the wall shear stress direction, i.e., 

w _ u, (5) 

U 0 U 

he noticed that the velocity profile can be represented by the two 
straightlines in the hodograph plane. This is known as the Johnston 
triangular model (Figure 1) . Variations of this model have been used by 
many. Eichelbrenner and Peube (2) considered that the cross flow component 
can be represented by a higher degree polynominal of the main flow com- 
ponent, and Shanebrook (19) used the fifth degree polynominal expression. 


« 
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Klinksiek and Pierce (8) produced a simultaneously laterally skewed 
boundary layer flow and showed validity of the polynominal expressions. 
Lakshminarayana, Jabbari and Yamaoka (9) used l/7th profile for the mainstream 
velocity and Mager's profile for radial velocity and predicted the momentum 
thickness and limiting streamline angle on a rotating single blade. The 
agreement between theory and experiment was found to be reasonably good. 

They also derived an expression for radial velocity profile in the nearly 
inviscid region. This is briefly mentioned in section 2.2 


2. 1 Equations of Motion : 

For the purpose of this analysis the helical channel (Figure 3), 
whose pitch is very small, is approximated by a flat circular channel 
with a leading and trailing edge. The calculation makes use of this 
approximation, but the asymmetry of the flow that would prevail in the 
case of a helical channel is taken into account (i.e., it is assumed 
that the trailing edge does not influence the leading edge flow). This 
model is similar to those of reference 9. 

The equations of motion in cylindrical coordinates (r, 0, z) 
rotating with an angular velocity Q about the z axis is given by. 


i iH + + 3w 4 . W 

r 89 97 8r r " 


( 6 ) 


u 9u 8u , 9u 

7 36 + v 37 + ” 77 


- — (2ftr - u) 
r 


JL 9p_ 1 9l u 

p r90 p 8z 


4 


(7) 
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u 3w , 3w , 3w 

“ -Tq + v -5— + W -5 

r 39 3z 3r 


(fir - u) 


1 ^ 1 3 t 

1_ _3p_ 1_ w 

p 3r p 3z 


Equations (7) and (8) are 0 momentum and r momentum equations respectively. 

u, v, and w are components of velocity along r, 0, and z directions. 

and x are shear s tresses in the direction of u and w. The flow is assumed 
w 

to be incompressible and steady relative to the blade. 


2. 2 Velocity Profile Analysis in the Outer Region : 

In the outer part of the boundary layer, the shear forces are considered 
to be negligibly small. For a fully developed flow (-^g- = 0) equations 7 and 
8 reduce to (assuming that the velocity in the external flow is fir) , 


3u , 3u w . 

v h w — (2fir - u) = 0 

3z 3r r 


3w . 3w u 

v -sr— + w -r — 

3z 3r r 


2fiu = — 
r 


If w is assumed to be a function of u, the following expression can 
be derived (9 , 21) for radial velocity 


w = c y 

U 1 M Ji w . Un 
U - u ) (C + -) 


The experimental data obtained for the boundary layer flow on a rotating 
helical blade are well represented by the relation 


U “ °* 3 / ^ (1 - 

U v u' 


as shown in Fig. 2, 


It is interesting to note that, while the 


velocity profile at the outer part of the stationary type boundary 
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layer flow is represented by the linear relation given by Johnston, the 
velocity profile at the outer part of the rotating type boundary layer 
flow is represented by the circular arc. 


2 . 3 Flow Near the Blade Surface : 

In a region very close to the blade surface, flow must be mainly a 
viscous one, since all velocities, including turbulent fluctuations, 
become zero at the surface. In this region, velocity components are so 
small that the inertia terms become negligibly small compared to Coriolis 
terms. In this region equations 7 and 8 become 


- 2Qw = V 


a 2 

0 u 


(13) 


2 2 

+ 2f2u = 1- V ” (14) 

r dz 

where v is the kinematic viscosity, U is the free stream velocity (ftr) . 

By eliminating the Coriolis forces from the above equations, following 
equation is obtained 


^ v + V (w + u - 0 (15) 

r dz dz 

It is considered that the cross flow component can be represented by a 
function of the velocity component in the direction of the streamline 


- = f (— ) 

u nr 


(16) 


Substituting this relation, equation (15) becomes 


9z 


r f + v + u + a?) 


Here, f' , f 1 ' represent the derivatives of f with respect to u/U. 
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It is generally accepted that the law of the wall does exist for 

three-dimensional boundary layer flows. Experimental verifications of 

the existence of the law were given by Pierce and Krommenhoek (14) 

for the stationary type boundary layer flow and by Lakshminarayana et al 

(9) for the rotating type boundary layer flow. It is considered that 

the velocity component in the direction of the streamline is represented 

as a function of z + = U z/v. 

T 


|p=F(z + ) (18) 

T 

In the region very close to the blade surface, it is considered that the 
flow is very little affected by the presence of the pressure gradient 
outside the boundary layer. Then, by the use of equation 18, after 
some calculation, equation 17 becomes 


— + ff ' 
U 


f" F' 2 = ■— (^J - ) 2 “ constant 


(19) 


r T 

The second term on the left-hand side is very small and hence let us assume 
it is zero for the moment. Later we will prove that this term is identically 
The first terra gives 


$+«• 


F" 


= constant C, 


= constant C r 


( 20 ) 

( 21 ) 


where and C ^ are constants, such that ^ 0* 

Equation (20) has a solution w/U = f = £ q (u/U)^ where £ q is a constant- (22) 
This also gives f" = 0 which shows second term on the left hand side is 
equal to zero. 

Therefore, substituting expression for (f) in Equation (19), 


zero . 
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F" 



(V 

V U ' 

T 


c 


2 


Upon integration twice, the expression for F is obtained. 


(23) 


F ~lF ' 2 ( 2+ / 2) 2 + C 3 z + + C 4 (2 A) 

Boundary conditions at the wall are 

dF + 

F - 0, —T-l at Z = 0 
dZ 

and 

F = Z + (1 - — 0 ■ ■ ~+ ( 77 - ) 2 Z + } = Z + (1 - BZ + ) (25) 

2(l+e o 2 ) r U T 

-2 -3 

where B is a constant with a value of order 10 - 10 for turbulent boundary 

layers depending upon the extent of cross flow. For small values of Z + , 
equation (25) approximates to 


F a Z + (26) 

This is identical to the expression given for the laminar sublayer of two 
dimensional boundary layer flows. 

The analysis developed above shows that the collateral region does 
exist in the rotating type three-dimensional boundary layer flows. The 
collateral region is defined as a region where the inertia terms are 
negligible and the viscous forces are balanced by the Coriolis forces. 

As the distance approaches zero, it becomes identical to the expression 
given for the laminar sublayer of two-dimensional boundary layer flows. 

In the present analysis, equation 11 is obtained in the outer region 
of the boundary layer where the shear stresses are considered negligible; 
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it is valid only in the region near the outer edge of the boundary layer. 
It is considered that majority of the boundary layer flow is represented 
by a region in which neither the Reynolds stresses nor the inertia terms 
may be neglected. Solution of this region can only be obtained by 
numerical integration of the entire equations of motion. 
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2. 4 Momentum Integral Solution : 

The method is an extension of Lakshminarayana et al's (9) method 
for predicting the boundary layer growth on a single rotating helical 
blade. The method described below includes the effect of the pressure 
gradients imposed by the external flow. The two momentum integral 
equations (one in the main stream and the other in cross flow direction) 
are reduced to two first order partial differential equations using a 
power law profile for the main stream velocity ^ Mager's profile for the 
cross flow and the skin friction relation is based on the boundary layer 
flow on a flat, stationary plate with modification to include the effect 
of rotation (9) . The imposed pressure gradients are obtained from the 
three-dimensional inviscid numerical solution of the flow field (15) . 

The resulting two first order partial differential equations are solved 
numerically by a finite-difference scheme. Details of the assumptions 
involved^ equations and the method of solution has been described in an 
earlier report (3D) . The two momentum equations in cylindrical coordinate 
systems are: 


a0 n 3U * 1 9U e 

~ST~ + IT vf" ( 26 12 - % > + (29 n 

e e 


30 12 * 2 C fB 

+ r IT + 20 12 + 6 2 ( k - D = — 


(27) 


9(r9 22 ) 

3r 



9U 

e 

3y 22 


1 3U 30 ,1 

4- O Q ^ — 4. 2_1_ 

* 21 U 39 90 


e u - (| - i) 


+ e 


'f0 


<h. 


"W 2 


^ o 


(28) 


*The equations developed and the solutions given in this as well as the 
next two sections have general validity; applicable equally to flow investi- 
gated here as well as the inducer flow investigated in Part 2 of this report 
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o 

where 8 , 8 „ etc. are momentum thickness, 6 , 6 are displacement 

1 J. A ^ 

thicknesses in main stream and cross flow directions defined as follows, 


'11 


/■6 

r (1 - dz 

o e e 


■x-I 


(1 - dz 
o e 


.6 


12 


/i u . w , 

- u" } if dz 

e e 


22 


- jV 

J o e 


) dz 


(29) 


f<5 

<““) dz 
o e 


21 


f 6 at 

“ Jou 2 


dz 


U 

and k = — , U is the velocity in the main stream direction at the edge 
fir 6 

of the boundary layer. 


The assumed main stream velocity profile is 

1 H-l H-l 

u .z.ti 

r • 1 ■ n 


(30) 


where H is the local shape factor and r) is nondimensional z coordinate. 


« 
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The assumed cross flow velocity profile is 


— = e (1 - n)' 

u w 


( 31 ) 


The skin friction relation used is given by. 


C 0 = 0.0582 R 0 1/5 (1 + 0.85 /el) ) (32) 


w 


where R 0 is the Reynolds Number based on relative tangential velocity and 
distance from the leading edge. It is proposed to use the follow 
Ludweig's and Tillmann's skin friction relation 


C ffl = 0.246 10" 0>678H R fl -°- 268 (l + 0.85 /E~Q) 

ID ^11 W 


(33) 


Use of the above assumed relations (equations 30-33) simplifies equations 
(27) and (28) in terms of three parameters 0^, and H. The resulting 
equations are , 


30 6 . 3U ' 3U a 

ir + < 2 + H > r 1 ar + 2Je w 0 n (1 + ^ 3r> + < e u e w > 

e e 


(34) 


3U 


+ ^ S w e il--2 

e 


C f0^ 


and 


3 < e w 0 ll> 

N - + 2NS 


6 11 3U e _ 2£ 9U e He 2 6 . + M (re 2 0 ) 

«r F +0 nT F w 11 9r w 11 

e e 


- H0„C4- 1) + 


11 v k 


w 2 


( 35 ) 
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where 

9 12 = n 2 (lln+7) = (30+14H) 

J e 6, . (2n+l) (3n+l) (3n+2) (H+2) (H+3) (H+5) 

w ll 

2n 2 (n+2) 16H 

(2n+l) (3n+l) (H-l) (H+3) (H+5) 

(36) 

6n^ = 24H 

(2n+l) (3n+2) (5n+2) (H-l) (H+2) (H+3) (H+4) 


= 9 3.1 n 2 = 2 

e w 0 11 (3n+2) (H-l) (H+2) 


L = 


e 6 

w 11 


e 


m = 


22 


e 2 e 

w 11 


We now have two equations and three unknowns (0^^ , and H) , since 

is assumed to be known and values of will be based on empirical 
relationships. Anyone of the following relationships for H ,can be used 
for the additional equation needed for the solution of the boundary layer 
integral momentum equations . 

(1) From the data of an isolated airfoil (VonDoenhoff and Tetrevin (26)), 


— + e 4 ' 68(H 2,975) {- r~ - -r— + 2.035 (H-l. 286)} = 0 (37) 

r dU r U do 

e f 0 


(2) Head's Enttainment Equation (27) 


a u a 

[U Q 011 H, x *) + f 0 U H 6 _^* - L ^ (U e e w e u ) = E(H 6 _ 6 *) (38) 


73F lu e °il 6-6 


where E(H^ ^*) is an empirical function obtained from two-dimensional 
experiments . 
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E(H 6 _ 6 *) = 0.0306 (H 6 _ 6 * - 3.0) J 
where H g _ g * = 2H/(H-1) 

(3) Integral form of Continuity Equation with external flow assumed to be 
two dimensional, inviscid and in radial equilibrium: 


(39) 


Iff <e n H) + r 





(r 


£ 

w 


6 11> - 



au 

dr~ E w 9 ll 


(40) 


as 

30 


(i + 


2a 

IT 


( i_k)) + s (1 _ 


2a ' l 
7T ; U 



2a 

TTU 


Iff <Y« 


e e 

where U ’ = fir - U and 6 = ^ - 0 , , , S is the semiwidth of channel and 

e e (H-l) 11’ 

a is a parameter to account for the velocity defect in the midchannel. 

* 

a - 0 (— ) if the two boundary layers from the adjacent blade surfaces have 
not merged. If the external flow has radial velocity component, as has 
been observed inside the blade passage near the tip when the two blade 
boundary layers and the annulus wall boundary layer have merged to form 
an "Interference region," a - 1. 

For small pressure gradient in the streamwise direction, H = 1.286 and 
equations 27 and 28 reduce to (for r ^ = 1.5 feet). 


H + 0.543 eR || + A (3.285TG + 0.543R ^2. + ^ (pM _ 1>5g8RG _ 1-805) 


- ~ ) - 0.0442 AR = 0 


(41) 


and 


3A 


2„ 3A 


ffff + °- 53 Y R w * A f<0 ' 3175 + 1.058RG) e " + e (2.4TC - + 1.06R 


dz 


w w 


3R 


ae 


+ + (1.073 - pp} + 0.0207 e w AR = 0 


(42) 
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where A = 8 n R 8 
He/ 

and R g is the Reynold 
the leading edge . 


11 


is the momentum thickness in the mainstream direction 
number based on peripheral velocity and distance from 


is the limiting streamline angle representing crossflow effect. 


1 

3U 

e 

TG = — 

30 

e 

p=5 1 r=> 
ii 

Q 

Pi 

3U 

e 

3R 

e 

U 

k = ” 
K ftr’ 

U : 
e 

ft is 

the 

> 

ii 

m ' 
+ 

0.85 

s 

R = r 

_ i 


represents a term due to radial pressure gradient. 


ft is the angular velocity of the rotor. 


w 


tip 


The boundary conditions are 


(1) At the leading edge A(0,R) = e (0 ,R) =0.0 

w 


(2) At A(0, 0.5) and £ (0, 0.5) are assumed to be known and are given 

w 

by the values for a single blade (9) 
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Effect of Laminar Region near the Leading Edge 

The flow near the leading edge is found to be laminar and its extent 

(0 ) has been approximately predicted by a critical Reynolds number defined 
c 

9 S 

as fHT 0 /v = 3 x 10 , where V is the kinematic viscosity of air. 
c 

Banks and Gadd’s (25) analysis is used to predict 6^ and £ y in the laminar 

region up to the point of transition. It is to be noted that this analysis 

is valid for. the laminar boundary layer over a single rotating blade, and 

so doesnot take into account the effect of imposed pressure gradients on the 

boundary layer. Since the actual pressure gradients (Fig. £) are very small, 

the neglection of this effect may be justified. For 0 > 0 , the flow is 

assumed to be turbulent and the turbulent flow equations (41) and (42) are 

used to predict 0,.. and £ . 

F 11 w 

2.4.1 Theoretical Results and Discussion 

The momentum integral technique described in the previous section has 

been used to predict the momentum thickness (0-^) limiting streamline 

angle at various r and 0 locations for (a) a single blade (b) four bladed 

flat plate inducer ( c) three bladed inducer with cambered blade. The first 

test case is used as a check on the previous analysis of reference 9. 

For the last case, the three-dimensional inviscid solution of Poncet 

and Lakshminarayana(15) is used to obtain the chordwise and radial variation 

of RG, TG, and K. The inviscid prediction of ip (static pressure rise 

s 

coefficient) obtained at mid passage of the three bladed inducer are 
plotted^ in Figure 5 and 6. The variation of K, RG, and TG are shown in 
Figures 7, 8 and 9 respectively. 
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The predicted variation of 0^ and for a single blade are shown 
in Figures 10 and 11 respectively. The momentum thickness predicted by 
this theory is in close agreement with those of reference 9. Use of 
new skin friction correlation seems to predict larger growth in momentum 
thickness towards the trailing edge (Fig. 10). The analysis developed 
here seems to predict the larger values of e towards the hub and the varia- 
tion with radius is not allowed for in the theory developed by 
Lakshminarayana et al (9). 

Figure 12 shows the chordwise (0) variation of the mainstream 
momentum thickness (0.^) for the three bladed inducer. It is seen that 
for 6 < 1.5 radians, the pressure gradients are small and the values are 
closer to those for a single blade, (y 1% of maximum value of 0 = 5 rad.). 

Further downstream (1.5 < 0 < 4 radians) 0 ^ is nearly constant 
and infact decreases at higher radii. This is unlike the single blade 
case where in 0^ continually increases with 0 reaching an asymptotic 
value at large values of 0. A study of equation (41) for 0^ shows that 
the effect of positive chordwise pressure gradients is to increase the 
boundary layer thickness, while the effect of positive radial pressure 
gradients is to decrease it compared to the single blade case. It can 
therefore be concluded that in the present case, the effect of radial 
pressure gradient is much larger than the chordwise pressure gradient. 

This is confirmed by large decrease in 0 near tip radius where the radial 
pressure gradients are found to be much larger than the chordwise pressure 
gradients . 

The experimental results of a four bladed flat plate inducer operated 
at 450 rpm and open throttle are also shown. in Fig. 12. They confirm the 
trend of the numerical solution for three bladed inducer but the 0 


2 0 


values are larger and fall in between the single blade case and the three- 
bladed case. This is due to the fact that the four bladed inducer has 
been operated at open throttle and therefore smaller pressure gradients 
compared to those of the three-bladed inducer. 

Figure 13 shows the chordwise (0) variation of the limiting streamline 
angle £ w at various radial locations. In the turbulent part of the 
boundary layer i.e., 0 > 0 c# £ w first decreases very rapidly with 0 up 
to 0 of the order of two radians then increases with increase in 0. 

At given chordwise (0) location, £^ decreases with radius (for values of 0 up to abou 

two radians) and increases with radius at larger values of 0. A comparison 

of £ w for three bladed inducer and single blade (Fig. 13) shows that 

the trend is similar in both cases except that there is no jump in £ w at 

0^2 radians for the single blade case. An examination of equation (42) 

indicates that £ is very sensitive to pressure gradients and its value 
w 

increases in the presence of positive radial and chordwise pressure gradients. 

As noted earlier both of these pressure gradients are positive for a three 
bladed inducer and hence confirm the earlier conclusion. 

The experimental results of £ w for the four bladed flat plate inducer 
lie in between the numerical results for the three-bladed inducer and 
the single blade values. The experimental results of £^ at the trailing 
edge for the three-bladed inducer are also shown in Figure 13. They fall 
in the same range as the numerical solution but the variation with 
radius is opposite. The numerical results show an increase in £^ with 
the radius, while the experimental results show an opposite trend. This 
clearly points out to the fact that the flow model (velocity profiles, 
shear stress etc) need to be improved to get accurate prediction of 



21 


£ and 0,,. However, it is quite evident from the study of the behavior 
w 11 

of the two boundary layer parameters 8 and £ that the presence of imposed 

1 -L w 

pressure gradients substantially change the boundary layer characteristics 
in a rotating channel. 


2.5 Asymptotic Solution : 

It is anticipated that far downstream in the blade channel, in the 
interference region (Figure 14), where the two blade boundary layers have 
merged an asymptotic solution should exist. Since the flow is assumed to be 
fully developed, all the velocity and shear stress components don't vary with 
0. Assuming the fluid to be incompressible and perfect gas, the flow as 
steady and applying boundary layer approximations, the equations of motion 
become _ 
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where U^, and are the velocity components in chordwise, radial and 
blade normal directions respectively, p is density, is the ratio of 
specific heats. From the energy equation, it is clear that pressure drop 
in 0 direction is due to shear gradients in the boundary layer and since 
for fully developed flow the latter don't vary with 0, the pressure drop 
is also constant in 0 direction. From the boundary layer approximation, 
the pressure variation in Z direction can be neglected. Assuming 
T = 3U e /9Z, the energy equation gives 3i/3z also independent of z. 

In other words, in the equation A3, the right hand side can be assumed to 
be a function of radius only. Differentiating equation 43 with respect to 
z and using continuity equation, we obtain. 


\N e 


y( F/O 
9a. ^ 


1 

A. 


J H Zh. 


cMe 

9a. 


9 tr* ) 


= 0 


(48) 


where F = fir - U . 

e 

Equations 46 and 48 can be solved by the method of separation of 
variables by writing. 


W e = A 1 (r) B 1 (z) 

Fr = (Or - U e )r = C^r) D^z) 

V » E (r) F (z) 
el 1 

The resulting equations are: 


dC C dA 
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(49) 
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Since 3t/3z varies in the radial direction only, 

2 2 

only. This gives d D^/dZ = constant a^ 


2 2 
3 U / 3z 
e 


varies 


radially 


(60) 


On integration, = a^(2S - z) z + a^. The velocity defect is maximum at the 
center of the channel and is given by 


U (r) = fir - - + 34.) a /sin cf) 

c r 5 T 


This gives 

U 

e 


„ . (S - (2S-z)Z) r~. r 

U + L — u a/sin <b 

c r 5 


(61) 


(62) 


(63) 


The other velocity components can be written in a convenient form as 

z + a. } (64) 

4 

2 2 (Z 3 -S 3 ) 

+ S(z Z - S Z ) - } (65) 

From equations 58, 64 and 65 it is seen that 


W = — • sin (j) { ( 2S - z) 
e r rv 


and 


V = C. c cos <f> {a.(z - S) 

e 1 j u 


W /fir - U = £ /sin <b where £ is a constant (66) 

e e e e 


and V /W r = (z-S) (l + 2/3 (z-S) 2 /(a 4- (2S-z)z)} cot <j> 

6 6 4 


(67) 


Equations (66) and (67) show that (i) radial inward flow is directly proportional 
to the velocity defect in the mainstream direction. 
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(ii) Blade to blade flow velocity varies directly as the product of 

radial velocity W , radius r and distance from the center of the 
e 

channel (z - S) . The predicted variation of U and W across the 

e e 

passage is in conformity with experimental observ tions in a 
three and four bladed inducer (Reference 22) . 


u~ 

2. 6 Turbulent Field Method 

In turbulent field method, the turbulence quantities are suitably 
modelled using either the turbulence stress equation or the turbulent 

f 

energy equation. This model equation along with the three momentum 
equations and continuity equations are solved simultaneously by finite 
difference technique. This method was proposed by Bradshaw (28) for 
two dimensional boundary layer. He later extended this to three dimensional 
(non-rotating) boundary layers (29) . Referring to Figure 16, the momentum 
and continuing equations referred to blade coordinates (x, z,r) are. 
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where C 

* 

P 

u,v,w 


curvature of the external streamlines, 
reduced pressure = x r) • (E5 x r). 

velocity components along the coordinate directions x, z and r 
(Fig. 16). 


T = 


3u 

3z 


3w 

3z 


3v , 2 

T = v 3 v' 

z 3z 


(72) 


^ x j^ z = components of fl in x and z directions respectively. 

Equations for the Reynolds stress in the streamwise (x) and radial (r) 
directions are (Ref. 27), respectively. 
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In the transport equations for Reynolds stress (equations 73 and 74), 
the effects of curvature and rotation are to introduce additional production 
terms(cu.+|£+ 2ft z ) ..w'v' in the mainstream (x) direction and (2CU.+ 2£^(-u'v') 
in the radial (r) direction. 

In addition curvature terms are introduced in the transport of stresses 
due to turbulence and viscous diffusion. Since the diffusion terms are of 
second order and curvature is assumed to be much larger compared to boundary 
layer thickness, its effect in these terms is neglected. Using Bradshaw's 
(28, 29) modelling f equations 73, and 74 can be expressed as: 
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(76) 


where T = maximum value of shear stress vector, 
max 


* The term cw (-v'w,' ) nn the right hand side of the equation is negligible 
compared to cu(-u'v')and hence is neglected. 
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where 


e = dissipation of turbulent kinetic energy 

(79) 


G x - (pV/p + 1/2 q 2 V r ) / | T max /P | 1/2 Wp 


where a , L , G , are Bradshaw's boundary layer parameters modified 

1 x x 

to include the effects of curvature and rotation. 

Equations 68, 69, 71, 75 and 76 constitute a set of five differential 
equations for five unknowns, u, v, w, X , T . These equations have to be 
solved simultaneously with proper boundary conditions. Bradshaw's numerical 
analysis of three-dimensional boundary layers is based on a functional 
relationship for a^, L, G based on two-dimensional boundary layer measure- 
ments. A preliminary attempt to modify these functions to allow for 
rotation and curvature is described below. 

A stress Richardson number, (R. ) to account for the effects of 

l 

x 

curvature and rotation, is defined as the ratio of production of a Reynolds 
stress due to curvature and rotation to the production due to mean shear. 
Richardson number in x direction is therefore 


(2fi + 
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3z v 
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(80) 
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and in r direction is 


R. = 
1 
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(2ft z + 2cu) 
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( 81 ) 


In most of the turbulent boundary layer flows, the production of a turbulence 
quantity is balanced by the dissipation and the convective and diffusion 
terms are in order of magnitude smaller. 


Production — Dissipation 


(82) 


For stationary boundary layers on flat surfaces we have 


Production due to Shear = Dissipation 
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For boundary layers with rotation, the production term in the shear stress 
equations are modified by (1 — R ) and (1 — R- r ) respectively. 
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Also the turbulent energy equation gives 
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Since dissipation is the same we have 


3 /r 

Tx„ 


Using equation (85), equation (88) can be expressed as 


. 0- R<xV 
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Equations (85) and (89) can be used to modify a., L and G as follows: 
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The Richardson number R. appears in equation (90) as the modification 
parameter. It can be related to a Rotation parameter R^ = e^S/r where 
is the limiting streamline angle at the wall and 6 is the boundary 
layer thickness as follows: 
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In rotating passages of axial flow turbomachines, ^ , c~ ”, 
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This completes the theoretical formulation. An attempt will be made in 
the future to solve these equations by finite difference technique. This 
method predicts not only the velocity profiles, but also the shear stress 
distribution inside the boundary layer. 
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3. EXPERIMENTAL PROGRAM 

3. 1 Facility Description : 

The blade was made of the transparent plexiglass of constant thickness 
of 0.5 inch. Both the leading and trailing edges of the blade were 
shaped to form a symmetric airfoil shape to assure the smooth flow at 
both edges. The plexiglass plate was heated in the oil tub, then put 
on the mold and was cooled to give an axial advance of 10 inches. The 
blade has a circumferential extent of 300 degrees. The details of the 
blade are shown in Figure 3. 

The reason for the use of the transparent plexiglass for the material 
is twofold. Since very little is known on the helical channel flow, 
application of various forms of flow visualization technique was con- 
sidered. The use of the transparent material as the blade would make 
application of the particle tracing method possible. The other factor 
concerns with the roughness factor the blade surface. The blade surface 
would be considered hydraulically smooth when the surface roughness has 
less than a certain value. The often used criterion for the hydraulically 
smooth surface is given by U^k/v < 2.5, where k is a height of the 
surface roughness. This condition gives the roughness to be less than 
0.02 in. which would be rather difficult to achieve when the blade was 
made of, e.g., fiberglass. 

Hub is essentially a circular cylinder, made of aluminum, of 18 in. in 
diameter and 23 in. in length, and 1/4 in. in thickness. Nose cone 
has a spherical shape of 18 in. in diameter, and was made of aluminum 
of 1.4 in. in thickness. First, the mold was constructed and molten 
aluminum was poured. This was later machined to give a smooth surface. 
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The rotor assembly was constructed by installing the four blades 
on the hub with an equal angular interval of 90 degrees. Each blade 
was fastened on the hub at eight positions by a screw of 3/16 in. in 
diameter. A 1/4 in. hole was drilled in the blade, and the screw was 
inserted from the hub. A circular hole of 3/4 in. in diameter was 
drilled on the blade surface at the end of the hole. This hole was 
filled by an aluminum disk of 3/4 in. in diameter and 1/2 in. in 
thickness, which was threaded. End of the screw inserted to the blade 
from the hub was in this aluminum disk, thus giving a firm support to 
the blade. 

Detailed drawing of the rotor assembly is given in Figure 4. In 
this Figure, the probe traversing unit is also shown. This is used 
for the measurement of the velocity distribution across the channel by 
traversing the probe while the rotor assembly is in motion. Seen in 
this drawing are, the probe holder which is seated inside the hub, 
location of the slots where the probe is traversed, the step motor 
which traverses the probe while the rotor assembly is in motion, and 
the slip ring unit which allow for the step motor to be controlled 
from the stationary system. 

3.2 Experimental Method and Instrumentation: 

3.2.1 Three Channel Pressure Transfer Device : 

In the past, pressure transfer device was installed inside the 
rotor assembly. Pressure measured in the rotating system was led to 
the pressure transfer device from one end and then to manometer through 
a stationary conduit which came from the other end of the pressure 
transfer device. The conduit was fastened to a stationary system to 
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keep it from rotating. This arrangement, however, had inherent dis- 
advantages. Due to vibration of the rotor assembly, excessive wear 
experienced by the device led to occasional damage to the unit. 

Since it was installed inside the rotor assembly, inspection of the 
device meant removal of the unit from the rotor assembly, which was a 
time consuming process. 

In view of these difficulties, a new pressure 
transfer devicewas constructed. Details of the newly constructed 
pressure transfer device are shown in Figure 4. The device has three 
channels, and each channel was made airtight by using the double sealed 
ball bearings. Pressure leakage along the cylinders was prevented by 
the use of the 0-ring and the plastic sealers. The three channel 
pressure transfer device was installed on a stand outside the rotor 
assembly. Pressure from the rotating system was transfered to the 
pressure transfer device through a shaft which was fastened to the nose 
cone, and its value was evaluated by the manometer reading. In order 
to reduce the possible interference from the incoming flow, the pressure 
transfer device was housed inside the streamlined cowling. 

In order to eliminate possible transmission of the rotor assembly 
vibration to the pressure transfer device, a coupling unit was con- 
structed (Fig. 4). The shaft attached to the nose cone was connected 
to the shaft of the pressure transfer device through a coupling which 
consisted of the conduit of 1.625 inches in length and two disks of 
synthetic rubber of .25 in. in thickness. Test showed that the vibration 
of the rotor assembly was not transferred to the pressure transfer by 
the use of this coupling unit. 
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3.2.2 Ammonia Transfer Device : 

Direction of the flow near the blade surface was determined by 
means of the ammonia streak method, i.e., by releasing a small 
amount of ammonia gas through the holes drilled on the blade surfaces. 
Ammonia transfer device is a means to transfer the ammonia gas from 
the stationary system to the rotating system. Details of this are 
given in Reference 5. The ammonia transfer device was mounted on the 
stand outside the rotor assembly. Details of the limiting streamline 
angle measurement are given in Section 3.6 

3.2.3 RPM Counter : 

The RPM counter consists of 4.5 in. diameter metallic disk fixed 
at the end of the rotating shaft. The disk has 60 notches at equal 
intervals on the circumference to allow light to pass from one side 
of the disk to a photocell on the other side. The protocell is 
connected to a digital counter which counts the pulses and thus the 
rmp to an accuracy of 1/10 rpm. 

Figure 17 shows a schematic diagram of the instrumentation used 
for the measurement of the velocities inside the rotor passage. Details 
of the hardware and technique are described in Section 3.7. 

All the measurements, with the exception of those described in 
Section 3.3, are carried out without the throttle at the exit. 

3. 3 Effect of Throttle : 

In order to investigate the overall performance of the four bladed 
flat plate inducer at various flow coefficients, a throttle was installed 
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at the exit. The throttle clearance was varied and the annulus wall 
static pressure, stagnation pressure distribution at the exit were 
measured. 

Figure 18 shows the static pressure distribution measured on the 
annulus wall with various throttle clearances. Numbers in the figure 
represent the throttle clearances in inches. The location of the 
leading edge of the blade is taken as the origin of the abscissa. It 
is seen that the annulus wall static pressure increases as the 
throttle clearance decreases. The negative pressure at the leading 
edge of the blade is considered to be caused by the blade blockage 
effect. Increase of the blade static pressure at the trailing edge is 
considered to be caused by the effect opposite to the blade blockage, 
i.e. , by the increase of flow area due to the absence of the blade 
thickness. It is seen that the static pressure is constant downstream 
of the rotor assembly. It is interesting to note that the pressure 
gradient on the annulus wall decreases as the throttle clearance decreases. 

The radial variation of stagnation pressure and angle distribution 
were measured 14 inches downstream of the trailing edge. The stagnation 
pressure rise coefficients at various throttle positions are plotted in 
Figure 19. Even at large throttle openings, appreciable pressure rise 
can be seen especially at the tip. This is caused by considerable 
boundary layer growth on the blades, and the shear pumping effect. The 
outlet angles, measured with respect to tangential direction, were 
derived using a. cobra probe. The flow reversal, near the hub, 
occurs at flow coefficient of 0.0592 or less (Figure 20). The performance 
of the flat plate inducer, plotted in the form of ij) - 4> curve is shown 


in Figure 21. 
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3.4 Blade Static Pressures : 

The flow in a channel differs from those around a single blade 
two aspects: (a) The pressure gradient that exists in a channel change 

the boundary layer behavior, local shear stress, limiting streamline 
angle etc. (b) The interaction between the boundary layers on the two 
adjoining surfaces gives rise to complicated profiles arising out of 
mixing. In view of this it is important to have a knowledge of the 
static pressure measurements on the blade. 

The blade static pressure measurements were carried, with no throttle, 
at 11 radii and at five tangential locations, 0 = 30, 75, 150°, 225° and 270 
on both leading and trailing surfaces of the blade. The results are 
plotted in Figures 22 and 23. 

In Figure 22, radial variations of the blade static pressure 
coefficient at various tangential locations are shown. It is seen 
that \p s takes nearly a constant value except in the region near the 
tip where it changes considerably, especially near the trailing edge. 

In Figure 23, chordwise variation of the blade static pressure 

coefficient at various radial positions are shown. It is interesting 

to note that the static pressure distributions on the leading surface 

behaves quite differently compared to that on the trailing surface. 

While lb takes the minimum values at the midchord location on the 
s 

trailing surface, it takes the maximum values on the leading surface. 

The existence of finite loading on the blades can be explained 
on the basis that the boundary layer on the leading and trailing surfaces 
grow differentially giving rise to asymmetrical flow around the blade 


chord. 
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3 . 5 Limiting Streamline Angle Measurement : 

The limiting streamline angles on a rotating body can be determined 
from the ammonia streak method described in Reference 5. Small amount 
of ammonia gas is released through the blade static pressure taps, while 
the rotor is in motion, and the trace of the gas is recorded on an 
ozalid paper. It is important to realize that, in order for the above 
procedure to be valid, the ammonia gas must remain in the collateral 
region, which means that the ammonia gas coming from the orifice on 
the blade surface must have a very low velocity. Typical traces of 
limiting streamlines derived by this method are given in Reference 5. 

Measurements of the limiting streamline angle were carried out 
at the predetermined tangential locations which coincided with the 
locations of the blade static pressure taps. At each tangential location, 
more than two measurements are made. Variation between the two measure- 
ments is found to be very small. 

In Figure 24 radial variations of the limiting streamline angle 
are shown for various tangential locations. It is seen from these 
figures that, except near the trailing edge, the limiting streamline 
angle on the leading blade surface behaves quite differently from that 
on the trailing blade surface. It is interesting to note that at 
0 = 30 degrees, the limiting streamline angles on the leading surface 
is nearly twice the value on the trailing surf ace . As 0 increases, 
this trend reverses itself, and at 0 = 150 degrees, the limiting stream- 
line on the leading surface becomes smaller than that on the trailing 
surface. Unlike the case for the rotating helical blade, the large values 
for the limiting streamline angles do not mean that the flow is laminar. 
The flow is found to be turbulent throughout the channel. 
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In Figure 25, tangential variations of the limiting streamline 
angle are shown for various radial locations. At r = 18 in., the 
limiting streamline angle increases as 6 increases on the leading 
surface while it takes a nearly constant value on the trailing surf ace 
For the radial distance less than 16 inches, the behavior of the 
limiting streamline angle is different from that of r = 18". As 
0 increases , it increases to reach a maximum value on 

the trailing surface while it decreases to reach a minimum value and then 
increases on the leading surface. This behavior is quite different 
from the behavior of the limiting streamline angle for the rotating 
helical blade which showed that the limiting streamline angle increased 
monotonically with increase in 0. 

It must be pointed out here that for the rotating helical blade 

(Reference 5) the tangential pressure gradient was considered to be 

zero, and the experimental investigation showed this to be the case. 

For the present investigation, experiment has shown that there exist 

the pressure gradients, as reported in Section 3.4. It is interesting 

to note that the behavior of the limiting streamline angle is similar 

to the behavior of the blade static pressure. This can be seen when 

Figures 25a~c are compared with Figure 23. This may be interpretted 

as the indication that the limiting streamline angle is closely related 

to the blade static pressure distribution. Thus, in order to evaluate 
the limiting streamline angle it is necessary that the pressure gradients 

be known. Dependency of the limiting streamline angle on the radial 

pressure gradient was reported in reference 9. Tangential variation of 

the limiting streamline angle can be predicted when the tangential pressure 

gradient is given. 
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3.6 Determination of the Wall Shear Stress : 

One of the most important problems in the boundary layer flow 
investigation is the determination of the wall shear stress. In two- 
dimensional boundary layer flows, the wall shear stress has been 
determined by various direct and indirect methods. The frequently used 
direct methods include the Preston tube method, the heated element 
method (HEM) and the floating element method (FEM) . 

Because of the presence of the centrifugal force, vibration and 
rotation use of FEM and HEM methods are severly limited in its applica- 
tion for wall shear stress measurements on a rotating blade. 

In this paper, the wall shear stress was evaluated by means of 
the Preston tube method. Validity of its use for three-dimensional 
boundary layer flows is discussed below. 


3.6.1 Law of the Wall : 

For two-dimensional turbulent boundary layer flows, it has been 
well established that, near the wall there exists a region where the 
velocity is represented as a universal function of z + by 


The function f is represented by 

= A 4n 2 + + & 


(93) 


( 94 ) 


where A and B are constants. This is known as the law of the wall. 

When a total pitot tube (whose diameter is smaller than the extent 
of validity of the law of the wall) used, the wall shear stress can 
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be related to the dynamic pressure AP recorded by the pitot tube 
through the following relationship. 


% d L 
f X) z 


F C 


At 3 


( 95 ) 


where d is the diameter of the Preston tube. 

This shows that the wall shear stress can be determined by the 
local x dynamic pressure, once the function F is known. Determination 
of the wall shear stress has been carried out by many (3,9,12). 
Experimental verification of the validity of this approach to three- 
dimensional turbulent boundary layer flows was carried out by Pierce 
and Krommenhoek (14). Results obtained by the Preston tube 
method were compared with those obtained by the floating element 
method which gave the direct measurement of the shear force, and good 
agreement was reported. This result is of considerable significance 
since this clearly shows that the two-dimensional similarity law was 
valid in the collateral region of the three-dimensional boundary layer 
flow. Since the thickness of the collateral region is about a few 
percent of the boundary layer thickness, the use of the Preston tube 
method is not considered appropriate when the boundary layer thickness 
is very small. 

Work done by Pierce and Krommenhoek clearly shows that the law of 
the wall does exist for three-dimensional boundary layer flow. For 
the derivation of the collateral relation for the rotating type 
boundary layer flow, it is considered that the law of the wall is 
valid. The investigations reported in Ref. 9 confirm this. 
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3.6.2 Applicability to Three-Dimensional Flows : 

It has been pointed out that, in order for the Preston tube 
method to be adequately applied to three-dimensional boundary layer 
flows, the diameter of the tube must be smaller than, or approximately 
the same as the thickness of the collateral region. In this case, the 
flow is regarded essentially as two-dimensional one and the law of 
the wall for the two-dimensional flow will be directly applied. Thus, 
in order to discuss the applicability of the Preston tube method to 
three-dimensional boundary layer flows, it is necessary to determine 
approximate thickness of the collateral region. 

The extent of the collateral region, however, is not determined 
clearly. In the present paper, the collateral region was defined 
as a region where inertia terms are negligible and the viscous forces 
are balanced by the Coriolis forces. In this respect, the collateral 
region may be considered to correspond to the laminar sublayer of the 
two-dimensional boundary layer flows. Therefore, the extent of the 
collateral region may be roughly considered about a few percent of the 
total boundary layer thickness. This severely limits the application 
of the Preston tube method to many actual cases. 

In the present investigation, it is found out that the directions 
cf the maximum total pressure obtained by the Preston tube did not 
agree with the limiting streamline angles obtained by the ammonia 
streak method. It is noticed that the difference in angles increase 
with the increase in limiting streamline angle, i.e., the difference 
in angles determined by two methods is large when the skewness of the 
flow is large. This trend is observed in the entire experiment. Even 
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though the extent of the collateral region cannot be determined 
defintely, the existing difference in angles determined by two methods 
clearly indicates that the diameter of the tube is larger than the 
extent of the collateral region. 

It is known that the pitot tube is insenitive to the direction 
of the flow up to a certain degree of angle. This means that the 
full amount (within + 1 percent) of the total pressure is recorded 
by the pitot tube when it is placed in the flow in such a way that 
its direction deviates from the flow direction by a certain degree 
of angle. This angular deviation, y degrees, is a characteristic of 
the individual pitot tube and is known as the sensitivity of the probe. 

In the collateral region, the direction of the flow is given by 

the limiting streamline angle a = arctan £ Q . As the distance from 

the surface increases, the flow angle, measured from the tangential 

direction, decreases. a is defined as a flow angle at the distance 

d (diameter of the tube) away from the surface. The direction of the 

flow changes from a degrees to a Q degrees within the distance d. When 

the pitot tube is placed in the direction oi^ + y degrees, it reads 

the full values of the total pressure for the flow whose direction 

lies between a to a + 2y degrees. This means that, for the flow 

oo 

whose direction lies between a and a + 2y, as far as the total 

o o 

pressure is concerned, is identical to the flow without skewness, 
i.e., two-dimensional flow. Though the value of y was not evaluated 
precisely, it is considered to be large enough to regard the two-dimensional 
calibration scheme to be an accurate one. This is the basis for the 
use of the two-dimensional calibration shceme for the wall shear stress 
evaluation of three-dimensional boundary layer flows. 
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3.6.3 Experimental Procedure : 

For the evaluation of the local total pressure, two different 
sizes of the pitot tubes were used. One was a tube of .028 in. outside 
diameter with 0.6 diameter ratio, and the other .065 in. in outside 
diameter with 0.7 diameter ratio. Both are approximately 1.5 in. in 
length. 

Measurements of the total pressure were carried out at five 
tangential locations, 0 = 30, 75, 150, 225, and 270 degrees on both 
leading and trailing surfaces of the blade. At each tangential 
location measurements were carried out at four radial positions, 
r = 11, 13, 15 and 17 inches. The tip of the pitot tube was placed 
approximately 1/8 in. upstream of the blade static pressure hole. 

The tube was rigidly attached to the blade surface by a cellophane 
tape. The other end of the tube was connected to the plastic tube. 

This was also attached to the blade surface by the cellophane tape. 

The plastic tube was connected to the pressure transfer device through 
the hole on the hub. The total pressure was evaluated by the manometer. 
At each position, starting from the tangential direction, several 
straightlines were drawn, with a five degree interval, on the blade 
surface through the point where the tip of the tube was placed. The 
tube was placed on the surface along these straigh.tlines , and the 
total pressure measured. When the total pressure thus obtained was 
plotted against the angle measured from the tangential direction it is 
seen that the total pressure takes a maximum value. This maximum value 
was taken as the total pressure at this position. 
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In Figure 26, the total pressures recorded by the manometer are 
plotted against the angle measured from the tangential direction. 

Since no centrifugal force correction was applied, the maximum total 
pressure is represented by the minimum value in these figures. The 
values of the limiting streamline angle obtained by the ammonia 
streak method are shown by the vertical dotted lines in these figures. 

It is seen that the difference between the limiting streamline angle 
and the angle that gives the maximum total pressure increases as the 
limiting streamline angle increases. It is also seen that the difference 
in angles decreases near the blade tip. This is considered to be the 
result of the increase in the boundary layer thickness. 

Using the blade static pressure obtained previously, the local 
dynamic pressure was obtained. Following the calibration procedure 
given by Patel (12) , the local dynamic pressure was related to the 
wall shear stress. 

All measurements were carried out using both tubes, and the results 
were compared. They are shown in Figure 27. It is seen that the two 
probes give almost identical values. With an intention to investigate 
the possible effects of the size of the tube, a pitot tube of .125 in. 
outside diameter with 0.6 diameter ratio and 1.5 in. long was used 
at some selected locations for the evaluation of the wall shear stress. 
The results agreed very well with those obtained by two other tubes. 

It is concluded that the size of the tube used in the experiment has 
little effect On the determination of the wall shear stress. This 


conclusion is identical to the conclusion obtained for the two- 
dimensional turbulent boundary layer flows.. 
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In Figure 27, local skin friction coefficients are plotted as a 
function of the Reynolds number based on the distance measured from 
the leading edge of the blade. The skin friction coefficient is defined 
by 



where U represents the local potential velocity outside the boundary layer 
and is assumed to be equal to fir. Since the velocity outside the 
boundary layer is not known, and since there exists the pressure gradient, 
it is not appropriate to compare the above results with the known 
expression given for the flat plate. 

Detailed discussion of skin friction measurements and correlations 
will be dealt with in a later section. 
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3. 7 Velocity Profiles Inside the Blade Passage: 

3.7.1 Method of Measurement : 

The mean flow velocity profiles inside the rotor blade passage 
are obtained by measuring total pressure and radial flow angle with a 
three hole disk probe. Measuring grid is 5 x 5 x 40 in tangential, radial, 
and blade to blade directions respectively. Figures 4 and 17 show the 
schematic diagram of the experimental setup. The radial and axial 
positions of the probe are adjusted manually while the blade to blade 
traversing (tangential direction) is done by a traverse unit installed 
inside the rotating hub. The traverse indexing unit is stationary and 
keeps the traverse wheel and the step-motor of the traverse unit locked 
except when indexed. The flow is traversed at five axial stations (Fig. 29) 
and five radial stations and 40 stations across the channel in the tangential 
direction. 

The disk type of probe has a number of advantages over a conventional 
yaw probe for measurements close to blade walls. It can be” easily aligned 
parallel to a blade surface, produces less flow distortion and has a wide 
linear flow angle measuring range (± 20°) as shown in the calibration curve 
Fig. 28. 

The mainstream velocity component u (parallel to blade surface) is 
obtained from measurements of the total pressure, flow angle and the static 
pressure interpolated linearly from blade static pressure measurements as 


• (Po~P s ) 



follows : 


Ut 
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( 96 ) 
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The radial component of velocity is given by 


W_ 

a 
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Displacement and momentum thicknesses are obtained by integrating the 
velocity component across the boundary layer using equations (29). Shape 
factor H is then calculated as the ratio of displacement to momentum 
thickness. 

3.7.2 Mainstream Velocity Profile : 

Figs. 30a to 30e show the mainstream velocity component across the 
passage at various radial and circumferential locations. From the profile 
distribution at the 1st slot (Figures 29 and 30a) i.e. 0 = 30°, it is 
evident that the flow has already become turbulent over most of the radial 
locations. The boundary layer is thicker on the suction side (as much as 
1.0 in. near tip). This may be due to inlet flow incidence because the inducer 
channel is being operated at open throttle. 

The boundary layer is thin and laminar near hub on the pressure side. 

(At hub, the Reynolds number based on peripheral velocity SLY and distance 
from leading edge is 0.8 x 10"*). 

The profile distribution at subsequent downstream measuring stations, 

(Fig. 30 b-e) shows that the boundary layer is turbulent at all the radial 
locations on both the surfaces and grows rapidly downstream. At measuring 
stations 4 and 5, (i.e. 0 ^ 210°) the two blade boundary layers have 

merged into each other near the tip. The mutual interference effect appear in 



49 


form of decrease in mainstream velocity component in the middle of the 
passage and the radial inward flow. The inviscid flow region disappears 
with the result that no definite boundary layer characteristics can be 
obtained in terms of the various types of boundary layer thickness 
parameters; ( b, & t On H e^c . ). It is obvious that in this region 

momentum integral technique can only give a qualitative picture of the flow. 
A modified viscid solution based on Cooper's finite difference scheme or 
Bradshaw's Method, modified and outlined in Sec. 2.6, is essential to 
predict this flow behavior. 

Comparison with Power Law Profile : 

In Figs. 31a-c the mainstream velocity component (u) normalized by 
the free stream velocity u^ is plotted against the distance from the 
blade wall normalized by the local boundary layer thickness (<5) . Comparing 

i/ n 

this with the Power Law profile of appropriate index ( 7 - - (»£/§) ) 

Ue v 

it appears that the experimental data fits this representation reasonably 
well. The agreement is poor near the tip because the boundary layer near 
the tip is affected by the tip clearance and the annulus wall effects. The 
The value of exponent varies from 1/6 to 1/11, except near hub at 0 = 90°, 
where the values are 1/4 and 1/5. Here the boundary layer appears to be 
in transition (Reynolds Number based on peripheral velocity and distance 
from leading edge 3 x 10^) . The value of exponent is different on suction 
side (from 1/8 to 1/11) than on the pressure side (from 1/6 to 1/8). From 
the plot of shape factor based on the exponent, it is seen that the shape 
factor increases slightly with 0 except near hub. The radial variation 
in the profile shape is large indicating that the radial pressure gradients 
are stronger than chordwise pressure gradients. This is true because the 
measurements were taken at open throttle, yielding mild chordwise pressure 


gradient . 
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Comparison with Cole’s Profile: In Figs. 32a-c, the mainstream 


velocity component u is plotted in the form of Cole's Law, 


XL _ L ( fttp t / 1 W IW) + C ( tj £*) 

U t " k ^ * 


(9fc) 


where u is friction velocity defined as U't = J , T 0 is 

the local wall shear stress, W is a wake function and Tf is a parameter 
which takes into account the effect of imposed pressure gradients. 

It is seen from the plots that there exist the law of the wall 
region ( ~ ■=. -L log -4- C ) up to Z. + ^ 2oo- 500 , 

depending upon the nature of the blade surface (suction or pressure) and 
the location of the measuring station. The mean value of constants in 
the law of wall vary in the range of 5.60 - 5.80 for the slope constant 
(1/k) and 4.5 to 5.5 for the constant C. The wake component W(y/g) of the 
Cole's Law appears to have different form than the sinusoidal form for 
stationary two dimensional flows and may be due to the fact that entrain- 
ment process which is dominant in the outer part of the boundary layer 
is different for the boundary layer on an enclosed rotating surface than 
for a stationary two dimensional surface. Value of 77 which is characteristic 
of the imposed pressure gradients is found to vary from -0.6 to 1.35. It 
is large and positive on the suction side and has a small negative value on 
the pressure side of the blade. It decreases with 6 up to 0 ^ 150° 
and later increases slowly. This is in accordance with the trend of the 
observed blade static pressure distribution, as shown in Figure 22. The 
value of 7T increases with radius continuously in accordance with the radial 


pressure gradients. 


51 


The experimental results, therefore, indicate that for the turbulent 
boundary layers in the rotating channel, the law of the wall fits the data 
well in the region close to the wall (i.e. inner layer 30 S i? < 300) 

and the power law fits the data well in the outerlayer. 

Distribution of Passage Averaged Mainstream Velocity : 

In Fig. 33 is shown the radial variation of mainstream velocity 
component averaged across the channel. The passage averaged value of u 
is found to be lower than the free stream velocity that exists in an inviscid 
flow (u = £Lr). The difference increases as the flow proceeds further 

O 

downstream of the channel, especially after 3rd station i.e. 0 > I 50 
The largest difference is near the tip where the two surface boundary 
layers have interacted and for 0 > 15 0* it extends up to midradius ( R~0 8) , 
This large reduction in mainstream velocity near the tip is probably 
due to the mutual interference of the boundary layers on the two surfaces 
of the annulus wall boundary layer, which are now completely submerged 
into each other and also due to the observed radial inward flow in this 
region. This is in confirmity with the increase in radial pressure 
gradients observed near the tip region. As explained earlier, the value 
of shear stress is not negligible near the center of the passage and the 
turbulence level in this region is maintained at the expense of mainstream 
velocity gradient in the normal direction. This results in smaller mean 
mainstream velocities. Also from the local continuity considerations (over 
a cross-section) since the radial velocity component is increasing down- 
stream of the passage, both in magnitude and extent, the mainstream 
velocity is bound to decrease as the flow proceeds downstream of the 
blade passage, till it leaves the rotor. 
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3.7.3 Radial Velocity Profile : 

Figures 34a-e show the observed radial velocity distribution across 
the blade passage. The radial velocity increases on approach to the wall, 
indicating the extent of boundary layer on the blade surface. The radial 
velocity should decrease very close to the wall, but since the measurements 
could not be taken for values of £+ 100®. (due to finite size of the probe 

and this traversing unit) , only the outer-region of the boundary layer 
flow has been scanned and the measurements are consistent with the expected 
behavior. 

Near the center line of the channel, the radial velocities are small 
of the order of (ft/sec). However, qualitatively, it can be said that 
the radial velocity in this region is higher hear the tip region where 
stronger radial pressure gradients exist. It is seen that for 0 = 30° 
and 90°, the radial velocity is positive and implies more negative radial 
pressure gradient than predicted by simplified radial equilibrium theory. 

(Fig. 35) 

The suction and pressure surface boundary layers near the tip have 

merged resulting in interaction with each other as well as the boundary 

layer on the annulus wall. The observed radial pressure gradient is more 

than that predicted by the simplified radial equilibrium theory 

( i.e. -Ms 2 / o_ U ) and this implies the presence 

SR ~ *R ( ^ Ut 

of radial inward flow near the mid passage. This is evident from the 
velocity profiles shown for 0 > 90° (Figures 34a-d) . The inward flow 
increases both in magnitude and extent as the blade trailing edge is approached. 


« 
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Comparison with Maker's Profile : 

The observed radial velocity component is plotted in Fig. 36a-c. 
as streamline angle (local radial to mainstream velocity ratio ) versus 
distance from the wall, normalized by the boundary layer thickness (£) . 

A comparison with theoretical Mager's Model based on the limiting 
stream line angle (obtained from momentum integral solution (Section 2.4) 
and Mager's Profile (equation 31) indicates that the experimental data 
fits this form quite well except near the hub or tip. The agreement would 
have been better if the observed limiting streamline angle data were used 
in Mager's profile. However, in the tip region where radial inward flows 
are present, Mager's representation fails because of the reversal of 
radial velocity (S type of profile) inside the boundary layer. In fact 
any form of similarity in velocity profile in this region is expected to 
fail, because of the interciation of annulus wall boundary layer and 
boundary, layers from the adjacent blade surfaces. 

Hodograph Representation : 

The radial velocity inside the boundary . Jayer is plotted in a hodo- 
graphic form (local radial velocity (W) against local mainstream component 
(u) both of ;which are normalized by the local free stream velocity at the 
edge of boundary layer) in Figures 37a-g. The raeasurmments at R = 0.56, 
0.71, 0.866 are shown in Figures 37a-d and for the tip radii (R = 0.93, 
0.986) in Figure 37e-g. It is interesting to note that in the region where 
the two boundary layers have not merged (R = 0.5 to 0.866) the measurements 
in the outer region fall in a narrow band, thus indicating the existence of 
similarity in velocity profiles. The discrepancies, especially in Figure 
3b, can be attributed partly to the experimental inaccuracy. The boundary 
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layers are very thin which thus introduces appreciable error in the 
measurement taken from a probe whose thickness is about 1/5 of the boundary 
layer thickness. Even near the trailing edge (Station if 5) on the trailing side 
the boundary layer near the hub is thin. 

The polar plot near the tip (37e-g) shows no similarity in 
velocity profiles. The interaction between blade and annulus wall 
boundary layers in this region makes the profile characteristics unique. 

Plotted in Figs. 38a and 38b are the hodograph plots of velocities 
on leading and trailing surfaces respectively. The data appears to fall 
in a very narrow band for the leading surface (except the data very near 
the hub or tip), and fits the parabolic curve given by, 


JbL = 0-3 
u 


U e ^ Uc ; 


This equation has been derived in Ref. 9 for a single rotating blade on the 
basis of negligible pressure gradients and shear stresses in the outerpart 
of the boundary layer. This shows that the pressure gradients on the 
leading surface in the mainstream direction are very small except near the 
hub and tip. The data for the trailing surface shown is in Figure 38b 
has more scatter. 

Polar plot of the velocity profiles near the tip surfaces are shown 
in Figure 39. The radial inward velocities are comparatively small at 
0= 210°. The polar plot of the outer layer seem to fall on vertical 
line. Since the' radial velocities in this region are small, no definite 
conclusion can be drawn with regard to the shape or trend, since the 
experimental inaccuracy may be large. However, at 0 = 270° the radial 
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inward velocities seem to decrease the mainstream velocity. The radial 
inward velocity subsist only in the outer layer adjacent to the midpassage. 

This phenomena in the interaction zone is different and intrinsically more 
complex than those occurring near the blade surfaces. In view of the 
limited data in the inward velocity region, no definite conclusion can be 
drawn. Predictions based on asymptotic analysis do fit the data well at 0 = 270°. 

In Fig. 40 a,b,c are shown a similar polar plot, the abscissa being the 
velocity defect in the mainstream. 

Passage Averaged Radial Velocity : 

In Fig. 41 is shown the radial distribution of root square mean radial 
velocity, averaged across the channel width. The plot shows that the 
radial velocity increases downstream of the passage, the increase is 
very large near the hub region. Near the tip region, the value is nearly 
the same at all tangential locations. Higher radial velocities in the 
boundary layer imply higher values of the limiting streamline angle and 
from the above trend we find that this angle should increase as the flow 
proceeds downstream and decrease with radius when the RMS radial velocity 
is normalized with respect to local mainstream velocity, averaged across 
the passage width. 

3.7.4 Boundary Laver Characteristic Parameters (©Ll, H, £„): 

In this section the boundary layer characteristic parameters 

(i.e. the momentum thickness in the mainstream direction (6^)> the shape 

parameter H and the limiting streamline angle (£■ )) derived from the 

w 

experimental measurements are compared with the predictions based on 
momentum integral equations (Sec. 2.4 of this report). 
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The observed values of 0^ (Fig. 42a, b) is found to increase slowly 

further downstream. Q ^ ncreases with increase in radius. The values 

on the suction side are higher than those on the pressure side. The 

predictions agree well with the measured values except near the leading 

edge, where the extent of laminar boundary layer growth is assumed to be 

known. The transition from laminar to turbulent boundary layer is assumed 

to occur instantaneously at a critical Reynolds Number 3 x 10"* (based on 

peripheral velocity and distance from the leading edge.) Compared to the 

boundary layer growth on a single blade rotating in stationary fluid, the 

boundary layer growth is much smaller in a channel due to the presence of 

radial and chordwise pressure gradients. The difference between the 

observed and predicted boundary layer growth is partly due to an error in 

estimating the extent of initial laminar boundary layer and partly due to 

an error in estimating the momentum thickness 0^ from measurements. 

The observed values of the limiting streamline angle £ (Fig. 42 ki & c) 

w 

first increase in chordwise (Q) direction reaching an asymptotic value at 
6^2 radians. 

The value of 6 decreases with increase in radius and reaches a 
w 

minimum at the tip. The magnitude of 0 is higher on the leading side 

w 

indicating that the boundary layer is more skewed and that the cross flow 

is larger compared to that on the trailing side of the channel. The 

predictions agree reasonably well except near transition. The value of 

in the laminar region is obtained from the approximate solution of Banks 

and Gadd and asymptotic value obtained by Cochran. The numerical results 

are liable to some error because the transition is assumed to be sudden 

and there is an order of magnitude of difference in the value of limiting 

streamline angle (£ ) for laminar and turbulent boundary layers respectively. 

w 
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An important difference between the single blade results of reference 9 
and the rotating channel is that the values are higher in a rotating 
channel and there is a radial variation due to presence of varying radial 
pressure gradient. 

The distribution of shape factor H obtained from the measurements 
is plotted in Figure 43. The value of H is higher on the pressure side 
(1.30 to 1.50) compared to the suction side of the channel (1.10 to 1.30). 
Its value is nearly constant near the midradius but changes rapidly at 
hub and tip due to the interference effects of hub or annulus wall boundary 
layers. The value is more at higher radii due to larger radial pressure 
gradient. A plot of shape factor based on power law fitting of the 
experimental data (Fig. 44) also shows that the shape factor is almost 
constant, except near hub or tip regions. Since the experimental data has 
been obtained by running the four-bladed channel at open throttle and 
hence with mild pressure gradients, the above results show that it is valid 
to take H constant in O direction in predicting 0 and in this case 
but it would be necessary to include its (H) variation for future experi- 
ments with stronger pressure gradients as outlined in section 2.4 of this 
report. 

3.7.5 Skin Friction Correlation : 

Skin friction coefficients derived from preston tube measurements 
are plotted against Reynolds number (based on mainstream velocity and 
local radius) in. Figure 45. In the earlier plot (Figure 27) the Reynolds 
number is based on the distance from the leading edge. 
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The skin friction coefficient for the rotating channel is shown 
compared with those of a stationary channel (fully developed) with the 
same mean velocity in Figure 45. The correlation for a stationary 
channel is given by 



C99) 


where dh = hydraulic diameter, which is equal to twice the normal 
distance between the blades at any radius. The skin friction coefficient 
for a rotating channel is higher than those of a stationary channel on 
both the pressure and suction surfaces. 

The skin friction coefficients derived from the velocity profile, 
using clauser chart, are shown in Figures 46 a and b. These values are 
in close agreement with the values derived from direct measurement 
(Figures 45 a and b). 
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4. FUTURE MEASUREMENTS 

Measurements of static pressure, limiting streamline angle and 
skin friction at the blade surface and the mean velocity profiles inside 
the blade passage will be performed on the four-bladed flat plate inducer 
at 450 rpm and at a throttle opening of (0.75 inch) corresponding to a 
flow coefficient ^=.05. In addition, hot wire measurements would be 
made to derive the turbulent intensity and turbulent stress components 
inside the rotor passage. The method of measuring static pressure, limit- 
ing streamline angle and mean velocity profiles has been described in 
Section 3.7. The skin friction would be measured by using hot film 
sensors on the surface at selected locations and the hot wire measure- 
ments would be done by using a three wire sensor described below. 

4.1 Hot Wire Measurements : 

Since the flow inside the inducer is highly three-dimensional, a 
three sensor hot wire probe will be used for the measurement of 
turbulent quantities. Also, the asymptotic analysis developed in 
Section 2.5 is valid if shear stress gradient varies radially only in the 
interference region. Since this region is far away from the solid 
boundaries, only turbulent shear stresses are expected to exist and 
therefore measurement of the turbulent stresses inside the rotor are 
necessary to check the validity of the above assumption. 

A special probe configuration shown in Figure 47 would be used 

to obtain the turbulent quantities from the direct measurement of 

the three mean sensor output voltages E^, E^ , E^ and the six products 

2 2 2 

of fluctuating voltages , e 2> e 3 » e 1 e 2 ’ e ]_ e 3 ’ e 2 e 3 ^ • ^ orienta- 

tion of the sensors is shown in Figure 47. x , z, r are the directions 
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of blade chord, blade normal and radial directions respectively. Sensors 
1 and 2 are in the xz plane and are orthogonal to each other. Sensor 1 
makes an angle with the x axis. Sensor 3 is located midway between 
1 and 2 and is at angle y^ to the xz plane (or plane of the sensors 1 
and 2) . 

The constitutive equations of the hot wire and a method of deriving 
the turbulent quantities from the hot wire voltage measurements are 
given below. 

Effective cooling velocity on hot wire 1 (shown in Figure 47) is, 

V e ^ a [(U-u-') SiN*,-tV+>>) (^so^h-K^O+u) Cos<X, + ( V+v) S!i ^'J Z (ioo) 


where k is a function of wire parameters i.e., £/d ratio. 



Since 

2 

k « 1 and U > u, u, w. 

V e 2 ~ 

t J SiNcrf, 

p 2 . (. u’smot ( -\y CoSo<, - 

VCj 

U SlNc<| 

V e| oi 

U SlN°(| 



( 101 ) 


( 102 ) 


Using Reynolds averaging; V = V + v where V and v are mean 

e l 6 1 e l e l e l 


and fluctuating components , 


Ve, , u Sin*, [ l+ £ CotV, - ^ Cot*, + | ] 


(103) 


U> S u ' Sin*, - -v' Cos 4 - \n 

c l 


w 


0 SiNof, 


(104) 
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Similarly we have, 



Ve t - u Cos^, [ 

« £ + ^ * 2 ( 0 asi? J 

(105) 

and = uCosoi, 

TJ SiN^, VL_ _ . 

yCosi, 


(106) 

l\nd 




Vc, =. ♦ w. C*r, + i l o 

Ccrt-ar.f- ] 

(107) 

and 




V e3 m U Cos 

r ( - w' siN 


(108) 

The hot wires 

to be used have a_ • = y^ = 45°, 

therefore , 


vi = n. r 

ft L 

z. a u J 



Vez t= 0 r 

■Tz L 

M-k*- - 2 + (5) 2 1 

1 I o 0 J 


(109) 

Ve, = £ [ 

i + bh + ^ -+-1 0= )*" 1 

2 0 2 u J 




From King's law for heat transfer in hot wire, the voltage E is given by 

Z "Z. n • 

Ej = E 0 i + Bj Vj ( i=.i y 2.ov3) (no) 


where n. and B. are wire constants, with n. = 0.45 to 0.50 and V. is 

ii l l 

the effective cooling velocity. Applying Reynolds' averages. 


Et = Eoi * b, v," 


win 


. . n. Bi V, 


2 Ei V,- 


Hi 


e,- ^ 


■Oi = 


— i — a. , 


hi ( E. i - Eoi ) 

2- E-, Vj 


( 112 ) 
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Substitution of equations 104, 106 and 108 in 112 results in, 


£, - 
K, 


- 2 . _ 




u'- v' +w' 

u 


uSv'-t W ‘ Z*l 

0 


( 113 ) 


£5 

k* 


. / 
a -vJ 


where 


K v 


- hi B,- 


2J5 E, V, 0 -" 0 


aw 


Solving for fluctuating velocity, components u, v', and w: 


u' = i C , e t + Aoc . e s ) 

ZCl+«0 K| Ki. k s 


o' _ j ( £ _ ^ 

Z K*. K, 


-J - c 

2.C»+ ot ^ kl 


(115) 




where * Z VJ j \j 


Therefore, various correlations u^ , v^ , w^ , uV, vw', uw can be derived 
from equation (115,1 A technique of deriving the turbulence intensities 
and correlations from the hot wire data is described below: 

From u’, v' and w^ in equation 115, making desired products and averaging, 

U 72 - a 1 r e> + §: + 4 o4 l ?L + a e£_ . 40c ( ^ . clc» ) 1 

4 f i +0 /Z) v 7. u -l ~Z * ■+ *+ \ — + - — J I 

^ J L Kl K »- k,Kx. k,K 5 V^Kj J 


(116) 


V*- 7 


f + ^ _ 2 e,e L I 
L k»*- K,K t J 


(117) 
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w' 2 - = 


J_ r eh 
4m*) x L k, 2 - 




z e > e *- _4 ( e i e s 

K i k, k x 


k t k s J 


UV 


4o^ L 



2 °c ^fk — z»<. e i g 3, 

k *-*3 k.,**, 



( 118 ) 


(119) 


vw' - — — r £ -2. ^£3 + z e » e * - 1 

4ct4o0 L k- ktKj < lS k, x J 


(120) 


u'w' 


J 


L 


e* 

*.£• 


^ e>e h - 20-^ (. -k ^ — 4<< e 3 "I 

* v Ka. k «*S k.j.VCj k ».J 

3 ( 121 ) 


Since all the quantities on the right-hand side of equations 116-121 are 

2 2~ 2 

known combinations of the measured values e^ , , e^ , e^e^, e 2 e 3 » e i e 3» 

the turbulent quantities , uv', vw and uw can be determined. A 

schematic of the circuitry consisting of a four channel constant temperature 

hot wire annemometer. Model 420A wide-band multiplier and available 

averaging circuits to be used in this program, is shown in Figure 48. The 

hotwire probe would traverse by the traversing mechanism described in 

Section 3 the signals would be transmitted to stationary system using 

MSD #6, a six channel mercury slipring unit, as shown in Data Transmission 

Assembly (Figure 4 ) . 
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4.2 Wall Shear Stress Measurements : 

Skin friction measurements at selected locations (5 radial x 5 
tangential locations) would be obtained by mounting a disk probe with 
two heated skin friction gauges as described in an earlier report 
(Reference 21). Measurements would be transmitted to a four-channel 
constant hot wire annemometer through the mercury slip ring unit. 

The probe will be calibrated in a wind tunnel with the use of a Preston 
tube to provide a known value for the skin friction stress. 
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FIG3: 4 BLADED INDUCER CONFIGURATION 


71 



u 

C 

B 

Q) 

>-< 

3 

co 

CO 

e 

O) 

X> 

o 

u 

B. 

to 

n 


CO 

•u 

o 

u 

u 

o 

< 4-1 

>, 

■u 


o 

CO 

< 4-4 


<U 

CO 

u 

0 ) 

> 

CO 

4-1 

to 

c 

•H 

4-1 

CO 

•H 



p44 






FIG6 : RADIAL A CHORDWSE. DISTRIBUTION OF BLADE STATIC PRESSURE 
ON 3 BLADE D INDUCER ( REF. 15) 




a n 


os 













KA Dif\N5 


















Radons 
















/N R.ADH- \N5 ALONG CHORD. 

























KEUFFEL a KSSER Cp. 






























FIG .12: TAN 



to X to TO THE CENTIMETER 46 1513 






























lO X 10 TO THE CENTIMETER 46 1513 



































die 


STATION A ftv CHANNEL 


asymptotic SOLUTION 

PoK- KOTOR. CHANNEL 


i distance normal, to WALL 


0-6 07 0-8 R=R C 0-9 

Radius r 


Fig 1 5 : Mature of Velocity Profiles in asymptotic region 



FIG 16 : COORDINATE SYSTEM FOR TURBULENT FIELD METHOD 
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fig IB: Annulus Wall Static Pressure Distribution 
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Nomenclature 


AX,N1,N2 

{ a >b , c , d , e , f g 

B 

C f 

E 

E 


Orthogonal coordinate system along hotwire axis 

Coefficient constants in derivation of hotwire 
equations (Appendix 3) 

Slope of hotwire calibration curve 

_2 

Skin friction coefficient (= T • 1/2 pw ) 

Output D.C. voltage of hotwire anemometer 

D.C. voltage of hotwire anemometer at zero velocity 


e 

F 

8 o 

h 

k 

L 

N 

n 

0(e) 

P 

R 

R.e.z 

R 

e 

r 

U 

V 

w 

u,v,w 

Z 

z 


= Output A.C. voltage of hotwire anemometer 

= Body force including viscous terms 

= Gravitational acceleration 

= Static head 

= Correction factor for the deviation from the cosine law 
= Constant used to relate hotwire voltage and velocity 
= Number of blades 

= Unit normal to blade surface 

= Terms of small order 

= Static pressure 

= Non-dimensionalized radius (= r/x^) 

= Rotating cylindrical coordinate system 
= Reynolds number (= Ur/v) 

= Local radius 

= Local blade speed (- Rtt) 

’= Resultant absolute velocity 

= Resultant relative velocity 

= Relative velocity components in the R,0,Z directions 

= Non-dimensional axial location (= z/r^) 

= Local axial location 
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Superscripts 


Subs cripts 
ax,nl ,n2 
eff 
m 

r, 6, z 
t 

1,2,3 


Nomenclature (con't) 

= Angle of the flow with inducer axis 
= Angle of the blade with inducer axis 
= Local tangential location 
= Kinematic viscosity 
= Fluid density 

= Normal Stress 

= Shear stress 

= Inlet flow coefficient 

2 

= Static head coefficient (= 2g h/U ) 

6 o t 

= Angular velocity of inducer 


= Time averaged or passage-averaged quantity 
= Vector quantity 
= Fluctuating quantity 


= Components along AX,N1,N2 directions 
= Effective cooling value sensed by hotwire anemometer 
= Refers to manometer values 
= Components along R,0,Z directions 
= Refers to values at the inducer tip 
= Refers to hotwire sensor values 
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1. INTRODUCTION 

The complexity of the flow in axial inducers makes meaningful 
predictions extremely difficult. The long and narrow flow passages 
between the blades generate thick three-dimensional boundary layers, 
often fully developed, with large turbulence levels. The secondary 
motions in these inducers are not confined to thin regions at the blades, 
but extend over the entire cross section of the flow region. 

However, in order to understand the extent of these three-dimensional 
effects, it is essential that there be a thorough understanding of the 
three-dimensional inviscid effects such as blade blockage, flow turning, 
etc. on the flow region. A method of obtaining the exact solution to the 
inducer flow field has been developed by Cooper and Bosch in Reference 1. 

This three-dimensional analysis employs an iterative numerical procedure 
to solve the equations of motion expressed in finite difference flow. 

The use of this procedure for the three-bladed Penn State inducer (Figure 
la) and the current analysis being undertaken to reduce the necessary 
computational time and incorporate the dominant viscous terms is reflected 
in Section 2. 

Section 3 is devoted to the experimental program being conducted 
in the three-bladed inducer. The experimental data is used to supplement 
the theoretical analysis being undertaken, and vice versa, in an attempt 
to gain a thorough knowledge of all aspects of the inducer flow field. An 
extensive survey of the blade static pressures is presented, as well as 
preliminary results on hotwire measurements within the rotating flow passage. 

A comprehensive report on these investigations is currently under 
preparation and will be submitted when completed. 


4 
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2. THEORETICAL ANALYSIS 

A thorough knowledge of all significant inviscid effects (blade 
blockage, flow turning, finite hub/tip ratio, etc.) and viscid effects 
(boundary layer growth, energy dissipation, etc.) is essential in the 
accurate prediction of the flow in all turbomachinery. Cooper and Bosch 
(Reference 1) have made an important contribution toward the three- 
dimensional inviscid flow solution. Application of this method of 
analysis to the Penn State inducer is given in Section 2.3. In addition, 
Section 2.3 includes modifications to the Cooper-Bos ch method which 
are attempted to help reduce convergence time of the solution and pro- 
vide a viscid solution capability based on empirically determined blade 
skin friction coefficients. A method of initializing blade flow parameters 
as input to the Cooper-Bosch method in a further attempt for a faster 
solution convergence rate is given in Section 2.2. 

2 . 1 Numerical Solutio n of the Inducer Flow Field : 

The availability of modern computers with large storage capacities 
and fast computation times greatly enhance the possibility of numerically 
solving the complete equations of motion. One of the early investigations 
in this area was made by Cooper and Bosch (Reference 1) for the case of the 
three-dimensional inviscid flow through axial flow inducers. 

The application of this method to the solution of the flow in the 
three-bladed Penn State inducer has been reported by Poncet and Lakshminarayana 
in Reference 2. Because of the iterative techniques employed in the 
Cooper-Bosch program, a large amount of computer time is usually required 
to converge to a satisfactory solution. An investigation of available 




mathematical methods to solve the four simultaneous nonlinear partial 
differential equations governing the inducer flow revealed that there 
was no alternate method which would solve the equations more efficiently 
or effectively with a minimum of programming effort that the exact method 
described by Reference 1. The next available approach for the speedier 
solution of the governing flow equations is the optimization of the input 
parameters of velocity and pressure which would allow faster convergence 
to the three dimensional solution. Cooper and Bosch have derived an 
approximate analysis solution to be used as an initial input to the exact 
program. This method derives the blade-to-blade average quantities using 
axisymmetric equations, then uses these quantities in a blade-to-blade 
solution of an integrated form of the scalar momentum equation in the 
tangential direction. The flow parameters derived by this method were 
used in Reference 2. However, a different method of developing the 
initial input flow parameters has been attempted in the present analysis 
(Section 2.2). 

In addition to the development of the flow parameters, a numerical 
procedure to generate the axial, radial, and tangential coordinates for 
the three-dimensional mesh grid points required in the Cooper-Bosch 
analysis has been developed, allowing for easier and more accurate 
flow geometry input to the exact program. 

The above features have been incorporated in a computer program, 
written in Fortran IV and run on an IBM 370/165, which produces the 
complete set of input parameters in the correct data format used by the 
Cooper-Bosch program. It is hoped that the present analysis, being 
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perfomied at Penn State, will provide a more accurate and less time con- 
suming method for initializing the necessary input flow parameters, thus 
significantly reducing the convergence time required in the solution of 
the three-dimensional exact problem. 

2. 2 Initial Input to the Exact Program : 

The input program developed in the present analysis is divided into 
two basic parts: the derivation of the thorough-flow mesh grid coordinate 

system used for the Cooper-Bos ch exact analysis, and the generation of 
initial velocity and static pressure parameters at each three-dimensional 
grid point describing the flow region. 

2.2.1 Grid Geometry 

The input data required for the generation of the three-dimensional 
grid coordinate system include: 

1. The blade profile coordinates of up to three reference 
blade sections defined at constant radii. 

2. The thorough flow geometry mesh coordinates (r/r^ and z/x^) 
in the axial-radial plane. 

3. The design values of air flow angles (Figure lb). 

Within the blade passage (between leading edge and trailing 

edge) , the values of the tangential coordinates are generated by inter- 
polation between the pressure and suction surfaces of the reference 
blade profiles. The tangential coordinates are determined at any radius 
(or, more precisely, the radii associated with the thorough-flow (R,Z) 
coordinates) by the interpolation between or extrapolation from two 
adjacent reference radii. 



-159- 


Tangential coordinates are expressed in radians, while axial and 
radial coordinates are nondimensionalized with respect to tip radius. 

2.2.2 Douglas -Neumann Analysis 

An initial estimation for the velocity and static pressure dis- 
tribution throughout the inducer flow passage is calculated by the two- 
dimensional Douglas-Neumann program described in detail in Reference 3. 

The technique employed by the Neumann program to solve a particular 
fluid flow problem is to use source distributions of appropriate strength 
on the surface of the blade profile in such a way that the flow normal 
to the surface of the body is either zero or prescribed. When the 
Neumann boundary condition is applied, an integral equation in source 
strength O is obtained 

- V ro • n = O(s) + | a(q) A(q,s) dq (1) 

'body 

where A(q,s) = n • V(q,s) and V ot is the onset flow. V(q,s) is the 
velocity at a surface point s due to a unit source at q. 

The solution for the general case of a lifting cascade at any angle 
of attack is calculated by superposition of three "basic flows" in 
such a way that the correct angle of attack is obtained and the Kutta 
condition is satisfied. The "basic" flows are: flow at zero angle of 

attack, flow at 90° angle of attack, and circulatory flow for each cascade. 
Superposition of solutions is possible because the potential equation is 
linear and the boundary condition on the cascade blade is homogeneous. 

In the Douglas-Neumann results, velocities and static pressure coefficients 
are normalized with the modulus of the average onset flow velocity 

V. . + V .. 

inlet exit 


2 


( 2 ) 
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2.2.3 Quasi -Three-Dimensional Modification 

The investigation of the flow around an isolated airfoil in a 
contracting or diverging stream is presented in Reference 4. This 
analysis provides a simple method of modifying the two-dimensional 
Douglas -Neumann flow solutions to account for the three-dimensional 
effect of the coverging or diverging streamlines. An expression for 
static pressure coefficient on the airfoil surface is derived as a 
function of channel slope, two-dimensional static pressure coefficient, 
and the Fourier coefficients of the blade profile. The analysis 
utilizes thin airfoil theory approximations and assumes that thickness 
effects are the same as in plane flow. The mean flow is assumed to be 
inviscid, steady, and incompressible, and the variation of channel height 
is assumed to vary linearly from leading to trailing edge while the 
length of the contracting section is assumed to be the same as the axial 
projection of the blade. 

In the present application, the expression for static pressure 
coefficient has been modified in an attempt to represent the flow about 
a row of two-dimensional infinite cascades. This quasi-three-dimensional 
approach has been applied to the two-dimensional results obtained from the 
Douglas-Neumann analysis of 2.2.2 for each of the three reference blade 
profiles used to describe the Penn State inducer. The effect of the 
converging channel as determined by the above analysis on the Neumann 
solution for the Penn State inducer is essentially to decrease the blade 
static pressure near the trailing edge. 
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2.2.4 Velocity and Pressure Input Parameters 

Once all values of static pressure coefficient are calculated at 
three constant radii reference blade profiles by the previous analysis, 
the values of ip s at the necessary grid locations specified by the three- 
dimensional Cooper-Bosch input geometry are interpolated between or 
extrapolated from these reference values. Thus, especially at the hub 
inlet where r/r^_ = .25 (the minimum reference radius is r/r^ = .50) } 
the calculated ip distributions are not exact and may appear distorted. 
One possible remedy is to run the Douglas-Neumann program for each of 
the flow streamlines used in the Cooper input, with blade profiles found 
by projection of the reference blade profiles to the required grid 
locations. This may be attempted at a later date. Nevertheless, using 
the present method of analysis, the static pressure coefficients are 
determined for each three-dimensional grid coordinate of the Cooper-Bosch 
system. 

Values of radial, axial, and tangential velocities are then derived 
with the aid of the Douglas-Neumann program, the quasi-three-dimensional 
results, and the grid geometry. Thus, all the flow parameters necessary 
for the Cooper-Bosch analysis can be calculated. The pressure and 
velocity values, along with other necessary input parameters, are placed 
by the input program onto 9 track, 1600 BPI (Bits per inch) computer 
tape for acceptance by the Cooper-Bosch three-dimensional exact solution 
program. 

Figures 2a, 2b, and 2c compare the pressure distribution derived 
in this analysis with design values, experimental results of Section 3.2, 
and previous solution results of Reference 2. 
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A comparison of the radial variation of axial and tangential 

velocities calculated by the above method with the experimental results 

of Reference 2 shows close agreement (Figure 2d). The agreement 

v 

between the measured ““and inviscid prediction may be fortuitious, 

U t 

since the axial velocity predicted at the same locations is considerably 
different from the measured values. 

A preliminary listing of the input computer program is given in 
Appendix 1. The listing of the Douglas Neumann program is given in 
Reference 3 and is not reproduced here. 


2. 3 Exact Three-Dimensional Solution : 

A method of obtaining the exact solution of the inducer flow field 
has been developed by Cooper and Bosch (Reference 1). This three-dimensional 
analysis employs an iterative numerical procedure to solve the equations 
of motion expressed in finite difference form. The nonlinear partial 
differential equations governing the flow in a rotating cylindrical 
coordinate system R, 6, Z are: 


0 Momentum: 
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Where w, v, u are velocities in axial, tangential and radial directions 
respectively. These equations are rearranged to give residuals which are 
reduced to zero by a relaxation procedure. 

In applying this method to the Penn State inducer, the flow is 
assumed to be incompressible, and a grid of 7 x 7 x 26 is chosen to 
represent the blade passage. The flow geometry is shown in Figure 3. 

The boundary condition to be satisfied on the hub, annulus walls, 
and the blade surfaces is W • n = 0 , where n is the direction normal 
to the channel boundaries. 

The first of the 26 axial stations corresponds to the upstream 
thorough-flow boundary where the initial conditions are applied. For 
the boundary value problem to be consistent, these initial upstream 
conditions must specify the three components of velocity and pressure, 
and the tangential velocity on the second axial station (which thus 
defines the swirl at the inlet of the inducer) . 

The last four axial stations correspond to the downstream flow 
through boundary, and extend about one-fifth of the chord length down- 
stream of the trailing edge. With W • n = 0 to be satisfied on these 
stagnation stream surfaces, the set of boundary conditions for the 
problem is complete. 

The Cooper-Bosch program has been suitably modified for use on 
the IBM 370/165 system at The Pennsylvania State University. For increased 
efficiency, the program has been compiled under a Fortran IV H level 
optimization procedure which shortens the time required for repetitive 
calculations, and production runs were submitted using the resulting 
object card deck. 
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2.3.1 Inviscid Solution 

Using the input parameters of velocity and pressure derived from 
the analysis of section 2.2, prelimianry results of the exact inviscid 
program confirm that a lower total rms (root mean square) residual is 
obtained than with the previous method of initializing the input 
variables. 

As an example, the final rms residual for the inviscid results 
of Reference 2 was 0.12450 after 68 relaxation cycles, whereas a 
similar value is obtained using the present analysis in 10 relaxation 
cycles. This amounts to a considerable saving in computer time. 
Twenty-five iteration cycles has reduced the rms residual to 0.10579, 
indicating that a faster convergence to the solution should be possible. 
Further investigation should be carried out to confirm the effective- 
ness of the input analysis as an alternative to the Cooper-Bosch approxi- 
mate solution method. 

In a further attempt to decrease the convergence time, the exit 
flow angle was allowed to change depending upon the tangential and 
axial velocities calculated at the inducer trailing edge. Since the 
exact downstream boundaries are not known in this type of problem, 
it was hoped that by allowing the downstream boundaries to adjuct 
themselves and thereby unload the blade trailing edge, a more exact 
definition of the downstream streamlines would result in lower rms 
residuals. Cooper and Bosch suggest a similar technique as a means of 
reducing rms residuals in their recommendations for future work. 

Since the extension of the stagnation stream surfaces downstream 
have been constructed to be uniformly periodic with a spacing of 2rr/N 
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(N being the number of blades) , the values of velocity and pressure at 
the downstream tangential channel boundaries would be equal. This con- 
dition is applied at the blade trailing edge after each iteration cycle, 

resulting in an outlet angle defined by 3 = tan ^ — . The change in 

w 

outlet air angle, if any, is then used to recalculate the downstream 
stagnation stream surfaces. 

This method also has the advantage of automatically forcing the 
Kutta-Joukowski condition for the blade pressure distribution to be 
satisfied. Changes made to the original Cooper-Bosch program can be 
seen in Appendix 2 and are concentrated in subroutine "Main". 

As yet, this analysis has been untested. 


2.3.2 Viscid Solution 

In addition to the attempts to improve the convergence of the 
exact inviscid solution, a method of incorporating viscid effects into 
the governing equations of motion has also been investigated. The viscid 
equations of motion are: 
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Z Momentum: 
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Where 
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Molecular viscosity terms have been neglected in these equations . 


Comparing these equations with the momentum equations used by 
Cooper and Bosch, the following expressions for FR, FT, and FZ (the 
exact program variables for viscous loss terms) can be given as: 


i 9t n 9t 
_ I r r0 rz 

FR “ " p r3(P + 


3z 


9 ° (S - O pq ) 

+ + — -El EE_] 


3r 


( 10 ) 


FT 


1 r ^ 69 + 3T 6z . 3T 9r 2 . 

p r30 3z 3r r 0r 


( 11 ) 



3x 


r38 


+ 


3a 

zz 

3z 


+ 


3x t 
rz rz . 

3r r J 


( 12 ) 


Since the stagger angle is very large, these viscous terms can be 
approximated at the blade surface by neglecting second order terms as 
well as normal shear stresses to the following expressions: 
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The distribution of wall shear stress is assumed to be linear across 

the flow passage from pressure surface to suction surface. The values 

of wall shear stress are assumed to be known from previous experimentation. 

Skin friction coefficient for a four blade flat plate helical channel 

is given in reference 5. The results, summarized in Figure 4, are 

considered to be valid for the three blade inducer under consideration. 

Interpolation of the curves in Figure 4 for a given blade surface grid 

1 -2 

location under consideration gives a value of shear stress x = C f • ^pW 

WR 

for the appropriate Reynolds number R g = — of the flow at that point, 
where W is the average relative velocity across the flow passage. Calcu- 
lation of T at each grid location within the flow passage allows the derivatives 
of equations (13)- (15) to be calculated by finite difference methods. 

The changes to the original Cooper and Bosch exact program necessitated 
by the inclusion of the viscous loss terms are made in subroutines "Main", 
"Dloss", and "itesid" and are shown in Appendix 2. 

The effect of the viscous loss terms in the exact three-dimensional 
analysis is, as yet, untested. 
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3. EXPERIMENTAL PROGRAM AND DISCUSSION 

The primary goal of performing the following experimental program 
is to investigate the flow characteristics such as velocity, turbulence, 
and static pressure of the relative flow inside a three bladed inducer. 

The importance of this experimental data for a better understanding and 
prediction of the flow is mentioned in section 1. 

A secondary goal is the determination of the effectiveness of triaxial 
rotating hotwire as a method of measuring mean and fluctuating velocities 
and turbulence stresses within the inducer blade passage. 

3. 1 Three Blade Inducer : 

The experimental investigation is performed on a three foot diameter 
axial flow inducer with three equally spaced blades. The test facility 
is pictured in Figure la. Design of the blades are by the mean streamline method 
of Wislicenus (Reference 6). The inducer is operated at 450 rpm, which 
is determined to an accuracy of 0.1 rpm by means of a photocell circuit 
with rotating calibrated disk and displayed on an electronic counter. 

Important parameters of the inducer are as follows: 


Number of Blades 

3 

Hub Tip Ratio at Outlet 

0.50 

Hub Tip Ratio at Inlet 

0.25 

Radial Clearance 

0.0625 

Inlet Flow Coefficient (Design) 

0.065 

Blade Chord at r/r^ = 1.0 

82.96" 

Blade Chord at r/r = 0.75 

63.18" 

Blade Chord at r/r = 0.50 

49.94" 
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Solidity at r/r = 1.0 

2.15 


Solidity at r/r =0.75 

2.21 


Solidity at r/r = 0.50 

2.61 


Reynolds Number Based on Tip Radius 

7.0 x 10 5 


Maximum Deviation of Camber Line 

r/r = 1.0 

.00637 

from Mean Streamline 

r/r = 0.75 

.010 75 

(— ) 

L max 

r/r^ = 0.50 

.020 

The design blade and flow angles are 

given in Figure 

lb. 


The use of the three bladed inducer for the continued experimental 
investigation defined in this report is a result of conclusions reached 
by prior analysis and detailed in Reference 2; namely, it has appreciably 
better performance than a similar four bladed inducer tested at the same 
flow coefficient. 

Blade static pressure measurements were obtained with use of hypodermic 
steel tubing of .063" ID imbedded in the blade at ten separate pressure 
and suction surface locations. The pressure measurements were carried out 
at five radii from taps of .063" diameter drilled at equally spaced 
intervals from tip to hub. The chordwise location of these stations 
are given in Table 1. The radial location of the stations are given in 
Table 2. 

Velocity and turbulence measurement stations within the blade passage 
have been previously used (Reference 2) and were constructed by cutting 
tangential slots in the hub wall at the locations described in Figure 3. 

The inducer is statically and dynamically balanced at facilities 
in the Garfield Thomas Water Tunnel of The Pennsylvania State University. 
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3. 2 Blade Static Pressure Measurement : 

An extensive experimental investigation of the blade static pressure 
distribution has been undertaken to help provide a check on prior theor- 
etical analyses and provides information useful in formulation of future 
theoretical development. 

3.2.1 Equipment Used : 

A schematic diagram of the blade static pressure test setup is shown 
in Figure 6. The equipment used to measure the blade static pressure 
distribution of the three bladed inducer described in detail in Section 
3.1 is as follows: 

1) Scanivalve - The scanivalve (Figure 7a) , a scanning type pressure 
sampling valve for measuring multiple pressures was mounted in the 
rotating hub section of the inducer. The scanivalve incorporates a fluid 
wafer switch for time sharing one pressure lead with up to twenty four (24) 

unknown pressures, and is stepped by a rachet-geared solenoid. 

2) Solenoid Controller - A solenoid controller uses push button 

pulse length feedback and increased drawing voltage to step the solenoid 
driven scanivalve. The controller was equipped with a 24-division 
indicator dial which allowed monitoring of the static pressure station 
under consideration. 

3) Three Channel Pressure Transfer Device (PTD) - A 3-channel 
pressure transfer device, constructed by the Penn State Aerospace 
Department, is used to transfer the static pressure measurements from 

the rotating reference frame of the three bladed inducer to the stationary 
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reference frame. Each channel is made airtight by the use of double 
sealed ball bearings, and pressure leakage is prevented by use of 0-rings 
and plastic sealers. The PTD is mounted on a stand outside the rotor 
assembly and is housed in a streamlined cowling to reduce any interference 
on the incoming flow. 

4) Slip Ring Unit - An eight-channel slip ring unit, constructed 
by the Penn State Aerospace Department, is used to conduct electrical 
signals from the stationary reference frame to the rotating reference 
frame of the scanivalve. Electrical continuity is provided by carbon 
brushes in contact with a rotating commutator aligned along the inducer's 
rotational axis. The slip ring unit is mounted on the pressure transfer 
device, and all electrical and pressure connections are transferred 
through a hollow shaft and flexible couplings to the nose cone of the 
inducer. 

5) Peripheral Equipment - A transistorized 30 volt D.C. power 
supply was used to provide voltage to the scanivalve and solenoid con- 
troller units. A micromanometer graduated in 0.001" divisions was used 
to measure the blade static pressure. 

3.2.2 Experimental Procedure and Techniques '. 

Flexible vinyl tubing of 0.063" inside diameter was used to connect 
the ten suction and ten pressure surface stations to the available tubes 
on the scanivalve (Figure 7b) . The vinyl tubing was also used to connect 
the collection tubulation of the scanivalve to the measurement channel of 
the pressure transfer device. The manometer was similarly joined to this 
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channel. Electrical connections from the scanivalve were transferred 


through the slipring unit to the solenoid controller and power supply 
(Figure 7c). Masking tape was used to cover the radial pressure taps 
not under consideration in a particular run. 

Once the inducer speed was adjusted to 450 rpm, the static pressures 
associated with the ten pressure and ten suction surface stations at a 
particular radial passage location could be measured by stepping the 
scanivalve through its range of operation. The dial on the solenoid 
controller would indicate which station pressure was being measured. 

Each step provided the blade static pressure of a different chordwise 
station. Pressure readings were measured on the micromanometer to 
0.001". Once all readings for a particular radial passage location were 
obtained, the inducer was brought to rest and the next radial passage 
location was considered by uncovering the static pressure taps associated 
with it and covering the previous tap locations. 

Since the blade static pressure measurements were taken on the rotat- 
ing blade, it is necessary to apply a centrifugal force correction to 


obtain the static head. If h is the height of the water column measured 

m 


by the manometer, the actual static head is given by 

u Pm u j. V 2 t 2 2, 

h = — h + — (r - r ) 
s p m 2g o 


where p is the density of the manometer liquid, r is the radius of 
m o 

the rotating shaft used in the pressure transfer device, and r is the radius 
of the static pressure tap under consideration. 
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3.2.3 Blade Static Pressure Results : 

Experimental results are plotted in Figures 8a-8e for the various 
radial passage locations. Comparison of these results with the exact 
analysis results of Reference 2 show close agreement, especially in the 
tip region (Figure 2a). The Cooper-Bosch program, however, predicts 
negative for the pressure surface near the leading edge for all 
passages considered, when in fact this result was' not noticed during 
experimentation. Near the blade tip, is positive near the suction 
surface leading edge and becomes negative at approximately 20 percent 
chord. The magnitudes of agree favorably, especially at the trailing 
edge where little discrepancy is noted. 

The values of \p obtained from the Douglas-Neumann approach of 
Section 2.2 agree favorably at the tip and mid passage, but is noticeably 
different at the hub, due perhaps to the method of extrapolation used 
in the construction of that curve, as mentioned in Section 2.2.4. The 
quasi-three-dimensional approach has been intended, however, only as 
a means of closely approximating input flow parameters for the Cooper-Bosch 
analysis and, thus, is shown in Figures 2a-2c only to indicate its 
relative proximity to the experimental and exact solutions. 

3. 3 Three Sensor Rotating Hotwire-Measurements : 

Initial feasibility investigation into the use of the hotwire 
anemometer in the rotating flow passage of the three blade Penn State 
inducer has been mentioned in Reference 5. One problem experienced has 
been the inability of the relatively crude slipring unit to maintain 
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the continuity of the hotwire circuitry. It did not allow suitable 
operation of the hotwire for extended periods of time to allow mean- 
ingful measurements to be obtained. The present study is an attempt 
to improve upon the accuracy and longevity of the initial investiga- 
tion and prove the suitability of hotwires in the measurement of the 
relative mean and fluctuating velocities in a rotating environment. 

3.3.1 Equipment Used : 

A schematic diagram of the rotating hotwire test setup is given 
in Figure 9. A detailed description of the equipment used in the measure 
ment of the relative mean and fluctuating velocities within the rotating 
passage of the three bladed inducer follows:. 

1) Triple Sensor Hotwire Probe - A subminiature triaxial probe 
designed for boundary layer flows was used in the experimentation 
(Figure 10). The wire is 3 ym diameter copper plated tungsten with 

a length /diameter ratio of approximately 300. The probe is attached 
to a specifically designed probe support for use in traversing the 
inducer flow passage (Figure 11a). 

2) Three-Channel Constant Temperature Hotwire Anemometer - 
Two dual channel constant temperature hotwire anemometers are being 

used to provide the three channel capability necessary for these measure- 
ments. The original circuitry is shown in Figure 12, but is modified 
slightly for higher output voltage levels. 

3) Mercury Slip Ring Unit - A ten-channel mercury slip ring unit 
was utilized in the rotating hotwire measurements conducted on the three 
bladed inducer. The slip ring unit exhibits the smallest and most stable 
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resistance in the transfer of measuring signals from the rotating 
electrical elements to the stationary electrical conductors. Contact 
between the rotating wires and the stationary contact screws is made 
through a round contact disc to which the rotating wire is connected, 
rotating in mercury. Triple distilled mercury is used to provide the 
greatest conductivity and the lowest noise level distortion possible. 

4) Peripheral Equipment - The A.C. voltage signals obtained 
from the hotwire anemometer are visualized on a four-channel storage 
Tektronix-600 oscilloscope. Instantaneous mean velocity D.C. voltage 
readings from the anemometers are displayed on a digital voltmeter. 

The fluctuating voltage (A.C.) signals are processed through 
a 5.0 KHZ low-pass filter driven by a 15-volt regulated power supply, which 
is used to cut off the high frequency noise which may have entered the 
circuit. 

Mean square values of the A.C. voltages are obtained by passing 
the signals through a true rms voltmeter and subsequently through a 
manually controlled signal integrator. The mean squared voltage is 
displayed on a digital voltmeter. 

A sum-and-dif ference circuit is utilized to obtain the sum and 
difference between the three hotwire signals needed for the turbulence 
intensity calculations (Section 3.3.2). 

A wave analyzer is used for sinewave signal generation to determine 
gains throughout the hotwire circuitry and the accuracy of the associated 
peripheral equipment. 

5) Calibration Equipment - A low turbulence calibration tunnel 
is used for the hotwire calibration. The horizontal wind tunnel has a 


4 
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test cross-section of 1 1/2" x 1 1/2" and operates within the range of 
air velocities of 0 to 300 feet per second (Figure lib). The calibration 
velocities were measured with a nonsliielded pitot tube and the micro- 
manometer of section 3.2.1. 

3.3.2 Experimental Procedures and Techniques : 

Measurements were taken at two axial stations, corresponding to 
approximately 33 percent and 90 percent of the blade chord (Figure 5). 
Various velocity measurements have previously been' performed at these 
stations (Reference 2) and thus a comparison of hotwire experimental 
results with these prior investigations are possible. 

Six radial stations (corresponding to r/r values of .973, .945, 
.890, .781., ..671, .548) at station 1 and five radial stations (r/r 
locations of .973, .945, .890, .781, .671) at station 2 are traversed 
at several tangential intervals in an attempt to get an accurate and 
detailed appraisal of the flow velocities, turbulence intensities 
and stresses in these regions. 

The three sensor probe is attached to a ten-inch section of . 
aluminum tubing and fixed in a particular radial and tangential locations 
by a coupling mounted in the inducer hub (Figure 11c). The probe 
is accurately aligned in the tangential direction with the aid of the 
guide vane attached to the probe's adjustable protection pin. Orienta- 
tion of the three individual hotwires was measured with respect to 
the (r,0,z) coordinate system by utilizing a linearly calibrated 
scale eyepiece in a 30-X Bosch and Lomb microscope. The direction 
cosines of this orientation was then calculated, as was direction cosines 
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of the two arbitrary normals to each wire. These values were used in 
the governing hotwire equations derived in Appendix 3. 

Tire experimental setup for the hotwire measurements is shown in 
Figure 13. With the probe ** position, the inducer is started and rotated 
to 450 rpm. The corresponding mean D.C. voltages of the three hotwire 
channels E^, are recorded, in addition to the A.C. fluctuation 

voltage values of e'^ 2 , e'^, e' ^ , (e^ + e ' 2 > 2 (e^ ~ e’ 2 ) 2 > (e ' ^ + e' 3 ) 2 , 
( e ' 2 + e ' 3 ^ » ( e ' 2 " e ' 3) 2 > (e^ “ e ' 3 ^ 2 * 

Mean voltages are obtained over an integration of 100 seconds. The 

inducer is then stopped, the probe is moved to another location, and 

the procedure repeated until the flow field is entirely surveyed. 

Station 1 and Station 2 are surveyed similarly, except that no turbulence 

measurements were recorded at station 1 . 

The resultant voltage measurements from the three channel rotating 

hotwire experiment are converted to mean velocities u, v, w and turbulence 
2 ~ 2 2 ~" 

quantities u' , v' , w' , u'v' , u'w’, v'w' from the appropriate cali- 
bration curves and the applicable equations derived from the hotwire 
orientation as explained in Appendix 3. 

The- data reduction was accomplished in a computer program written 
to solve, for all flow stations considered, the resulting three simultaneous 
mean velocity equations and six simultaneous turbulence velocity equations. 
The IBM 370/165 of the Penn State computation was used in this task. 

3,3.3 Mean Velocity Profiles : 

The mean velocity profiles of W/U.^ , u/U^_, v/U^_ and w/U^ are plotted 

versus percentage of passage width for station 1 in Figures 14a-14d 
and for station 2 in Figures 17a-17d. The variation of these profiles 
with various r/r t can readily be seen. 
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Measurements at Station 1 : 

Total Relative Velocity - Figure 14a shows the tangential variation 
of total velocity W/U . Comparison of this profile with the results 
of Reference 2 obtained by pressure probe measurements show good agree- 
ment at two similar radii (Figure 15). A velocity deficiency is noted 

in Figure 14a at approximately 55 percent 
passage width for all radial stations, but it is especially pronounced 
near the tip. This seems to indicate a concentration of high turbulence 
intensities in this region and, indeed, qualitative measurement of the 
A.C. fluctuating voltages confirm this fact. From the passage measurements 
of Figure 14a, the suction surface boundary layer can be discerned, appear- 
ing to grow in thickness near the tip. The effects of tip losses are 
apparent from Figure 16 where total relative velocity is plotted versus 
radius. The deviation of the flow from ideal conditions is apparent 
near the tip. As no evidence can be detected of the pressure surface 
boundary layer, it indicates that the suction surface boundary layer 
is larger than that of the pressure surface. It should be remarked 
here that no measurements are taken close to the wall. Since the blade 
element is not radial, the probe could not be located very close to the 
wall. 

Axial Velocity - Figure 14c shows axial velocity w/U t versus percent- 
age of passage width. The general trend for the tangential variation of 
axial velocity at the radii considered shows an increase toward the 
pressure surface. The radial variation of the axial velocity shows 
the largest values near the hub, decreasing consistently towards 
the tip. This tends to indicate the effect of blade blockage on the 
axial velocity distribution. 
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Raclial Velocity - Figure 14d depicts the tangential variation of 
the radial velocity u/U^,. Large values are shown to exist near the 
suction surface at radii close to the hub, indicating higher radially 
outward flow in this region. Near the tip regions, the profiles show 
a contour change at mid passage which may be interpreted as an indication 
of radially inward flow. 

Measurements at Station 2 : 

Total Relative Velocity - The tangential variation of total velocity 
W/u t is shown in Figure 17a. Again, as in station 1, a velocity deficiency 
is noticed at approximately 40 percent passage width. A growth in the 
dimensions of the eddy near the trailing edge can be seen by comparison 
with Figure 14a. The relative velocity distribution from R = 0.781 to 
0.973 are similar with very nearly same average values. Furthermore, 
the average values (w/u t ) are considerably less than the design values. 

Axial Velocity - Figure 17b gives the tangential variation of 
axial velocity w/U^. Overall magnitudes, of course, are higher than 
station 1 due to the converging annulus. Again, as in station 1, the 
radial distribution of axial velocity show larger values toward the 
hub to indicate the continuing presence of the blade blockage effect. 

In comparison, the opposite trend was found in Reference 2 at downstream 
locations where no blade blockage effects should be present. The pressure 
surface boundary layer is evident near the hub and displays a greater 
thickness in this region than near the tip. 

Radial Velocity - The variation of radial velocity u/U t across 
the passage is shown in Figure 17d. The overall magnitudes appear 
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large r than at station 1. The radial variation of velocity show higher 
values of radial velocity near the hub, with smaller magnitudes progress- 
ing toward the tip. No evidence of radially inward flow in noticed at 
midpassage, although the large change in velocity at approximately 25 
percent passage width near the tip indicates perhaps some change of 
direction may be occuring there. 

3.3.4 Turbulence Intensities and Reynolds Stresses : 

The turbulence intensities / , 2 /W, / = T2’/W , / , 2/W are plotted 

u v w 

in iso-contour form in Figures 18a-18c. The contours for all three 
intensity components are essentially the same, showing a pocket or 
core of high turbulence centered at 40 percent passage width and this 
agrees with the total velocity deficiency noted in Section 3.3.3 in 
this region. Turbulence levels are generally higher near the tip, 
whereas magnitudes of the hub are generally constant across the blade 
passage . 

Relative magnitudes of turbulence intensities rank /— r/W highest, 

i ^ 

w 

with / — r/W and / — r/W following in order. Maximum turbulence intensities 

f ^ | ^ 

U V 

of 25 percent are measured, which are considered high for the linear 
assumptions made in the derivation of the hotwire equations. Hence, 
some error is expected in the evaluation of intensity from the hotwire 
equations. 

2 2 

The isocontour Reynolds stress distributions, u'v'/W , u'w'/w 
2 

and v'w'/^ within the blade passage are given in Figures 19a-19c. 
Concentrations in stress intensities are similar to the previous turbulence 
intensity contours, but with slightly different passage locations. 
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4. CONCLUSIONS 

A brief summary of conclusions reached by the investigation recorded 
in this report are as follows: 

1) The input program derived in Section 2.2 appears to provide a 
better initialization of velocity and pressure parameters needed for the 
Cooper Bosch exact solution. Evidence of this is shown by initial exact 
solution inviscid results which indicate a reduction of total RMS residual 
after a similar number of relaxation cycles from results of Reference 2. 

2) Modifications to the Cooper Bosch program to automatically unload 
the trailing edge station and incorporate dominant viscid effects have been 
made in subroutines "Main", "dIoss", and "Resid". These modifications have 
yet to be tested, but significant .improvements are expected and will be 
presented in the comprehensive report to follow. 

3) A method to obtain velocity and turbulence measurements in the 
rotating blade channel has been tested. The triaxial hotwire probe in rota- 
tion utilized in this study have yielded satisfactory velocity profiles and 
turbulence intensity measurements. Comparison of velocities derived from 
rotating hotwire with those derived by rotating pressure probe measurements 
of Reference 2 show good agreement. It can be concluded that hotwire 
anemometry can be an extremely useful tool in the experimental determination 
of relative flow parameters in a rotating environment. 

4) Turbulence levels within the blade passage, indicated from the 
experimental results of this report, are generally high near the tip regions 
and, in particular, a growing core of high turbulence is evident at mid passage. 
This necessarily indicates the significant extent of flow mixing due to 
boundary layer interaction, radially outward flow, annulus wall effect, etc. 
which is prevalent within the long narrow passages of the inducer channel. 
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5) Overall agreement between all velocity components and pressure 
distributions obtained by present experimentation, prior experimentation, 
present theoretical analysis, and prior theoretical investigations is good. 
Detailed comparison will be presented in the later report. 
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APPENDIX 3: Derivation of Equations for A Three- 

Sensor Probe 

A3 . 1 Direction Cosine Method for Non-Or thogonal Probes : 

Consider one hotwire senor with respect to the (K.-e-i) coordinate 

system. 



Hotwire (1) has an orthogonal coordinate system associated with its 
position (A^\NL X N1). This ( /V^NL.NL ) coordinate system can be trans- 
formed to the ( ) coordinate system by: 

b, + c, ur 
W mi * CKi rr 4- ^ 

W/wl £ °^i ^ + b 3 va. C A <*r 


(1) 

( 2 ) 
(3) 
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where a^ , b^, c^ , etc. are the associated direction cosines between 

and ((3.-0-3L ). The effective cooling velocity is known 
(Reference 7) to be 

VvU +■ C 


or 


VO 




- (vw, * hx + vl 


(*) 


Considering instantaneous quantities, equationS(l) - (3) can be rewritten as 

V0**= +- a, [ir * ^ ,N ) + ^ c, (w-+ur} 

Vv) M »' 6 \N^ <Kx{*-+*\ 


(5) 


( fi ) 

VvJ^- Vi^ r “ + V | -v- ]o 3 \»L 4- W ^ *- c 3 (* + ^ 

; (7) 


Substituting equation (5) - (7), into (4), we get 

V\)^- \M«^+ r Iff’**) V \>* v^\ 


(8) 


— V \ 

+-Cq (u x + uJ X 4-Zw- +- <M (/***-* ajV + U. ror -r ^ J 

-V“ Q.^ (ax-ut Jr /N^ur + ^ \rT -V Jr ui W* -4- T*. W* + W*U.y^ 


where the constants are defined as 


M - Z (fc*C x v lc 

C M - +“ ^ 4 " ^ *fv^ * Z (b v C z 4* \>3^ + £ l, 


X / 2 6 

Equation (8) can be linearized by (1 + 6) = 1 + -j > where 6 is considered 

negligibly small, to the following expression: 


VO 




to the following expression: 

(gw) f < 



182c- 


/ 


/ 


/ 


Taking the time average of equation (9) , we get 

'TC, 

Subsequently, the difference between equation (9) and (10) is the 

i 

fluctuating velocity W ^ 

1 \i — 1 d-H * JLq * 

A3. 2 Relating Voltages and Velocities ; 

Apply King's Law for instantaneous velocities (Reference 8) 


then 


-v- 2. 4., *»- ^ - x-o 


v £= e 6 + e>'C?n 


W»4 \M^ 


,2 


VJ. 


Assume e and - iii r are small. Linearize the equation to get 

w > - W, 


i • & site, 


or 


w 




I 

Rewriting, Vd - L a , then 


i . «• i - 

\N, » L, ^ 


Vf -- t* <? _J 


AnO 


Vl^VO^. 5 

w[\m^ 


The value of fi ,- L 


**■» 1 




are obtained directly from the 


hotwire measurement of Section 3. 3; 2. The values of ^ Jif 4-^ 




are derived as follows 


(10) 


(ID 


( 12 ) 


(13) 


and 



I £2 d 


- (f 4] 

* T [C*>4 - K 7 ^-] 

^•i^a - 4 ^ ( £ l^- 5 -’.') - U^V)] 


Where the mean-squared voltage within the bracket are obtained directly from 
the sum and difference circuit measurements of section 3.3.2. 

A3 . 3 Mean Velocity Calculation 

Equation (10) for hotwire channel 1 and similar equations derived 
for channels 2 and 3 form a set of three linear equations in three unknowns. 


^I| — a 5“ + +* Cj 

V\)^ * Tr Jr b k t\. -v- C(„ ur 

V0 3 - Ol-, at +- ^ "^r 


where constants a,., b,., etc. reflect the appropriate combination of direction 
cosines, k factor, etc. 

Values of , V\J^ are known from application of equation 

(12) to the D.C. hotwire voltage obtained in Section 3.3.2 Solution of 
equation (15) thus give UL , A3" ^ for each measurement location 

considered. 


A3. 4 Fluctuating Velocities, Turbulence Intensities, and 
Turbulence Stresses 

Equation (11) is squared and time-averaged to give the following 
expression: 




u!at' 

+ 0 ( 6 ) 


12 2 e 


or 


V\f, =. a.o fT + L va! 4- c # *r‘ +■ JLo i\V 




.i 


t 


4- at'w-' + 


fj 


\a! oj* 


(16) 


v where a„, b Q , etc. are the appropriate constants, 
o o 

Similar equations exist for hotwire channels 2 and 3. Equation (11) 
can then be multiplied by the corresponding expression for V\J* and averaged 


to obtain 


(17) 


AT* + \j<j uj +■ ur*"*’ *r ^ (aIaT 

4- ^ at'vat' ulvr* 

where ^ ^ etc. are the appropriate constants. Similar expressions 


can be derived for and 


.i . » 


The values of wj , Vs)^ , Vsl,Vsl4, , ^i'Wi and 

are obtained from the analysis of Section A3. 2 

Thus the equations (16) and (17) form a set of six equations in six 
unknowns, which can be solved simultaneously to give the quantities UL* , (\J~ | 


U)’ 


U.aT, ( 




i. i 


AT VaJ and U. Ur . 
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